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Abstract 

We derive the Euler equations from quantum dynamics for a class of fermionic many-body 
systems. We make two types of assumptions. The first type are physical assumptions on the 
solution of the Euler equations for the given initial data. The second type are a number of 
reasonable conjectures on the statistical mechanics and dynamics of the Fermion Hamiltonian. 
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1 Main notations 



fj, Typical index labeling the five conserved quantities. 

j Typical index referring to time, j = 0, and space, j = 1,2,3, compo- 

nents. 

w w = (w^) = (wj x ),n = 0, ...4,j = 1,2,3, are the components of the 

current densities. \x = is the particle current, \x = 1,2,3, indexes the 
three components of the momentum current, p = 4 is the energy current. 
j = 1, 2, 3 refers to the three spatial directions of the current. 

W ... All quantities denoted by W*, are understood to be given by Z 3 w*, for 
any subsript *. 

w .... Underlined vectors have an additional component, referring to time, or, 
in the case of the currents, the conserved quantities. 

V .... V = (d t , V) is the four-component gradient, inclding the derivative with 
respect to time. 

u u = (u°, • • • , u 4 ) = w := wo are the five conserved quantities: particle 

number, three components of the momentum, and energy. 
.... u x = n x is the quantum observable for the particle number density at 

microscopic space point x. 
u x .... u x = (u x ,u x ,u x ), the quantum observables for the three components of 

the momentum density at x. 
u\ .... u\ = h x is the quantum observables for the energy density at x. 
A A = (Aj ) are the classical currents appearing the RHS of the Euler 

equations (/x = 0, ■ ■ ■ , 4, j = 1, 2, 3). They are defined by 

A°j = 

A) = dijP + qiqj/qo 
A] = \q A + P)l % • 



A .... A = (Aj) are the classical currents augmented with the conserved quan- 
tities (j = 0, vanishing current in the time direction): A = q = Ao- 

q q = • • • , q 4 ): q° = p is the classical particle density, q 1 , q 2 , q s are the 

three components of the classical momentum density, and q 4 = e is the 
classical energy density. These quantities are a function of macroscopic 
space and time. This notation is also used for the expectation value of 
these quantities in a quantum Gibbs state. 

q q = (g 1 , </ 2 , </ 3 ) are the three components of the classical momentum 

density. 

v v = (\/p is the classical, macroscopic mean velocity per particle. 

e e = e/p is the classical, macroscopic mean energy per particle. 



P(e, p) The thermodynamic pressure as a function of the energy and particle 
densities. Appears in the Euler equations and is defined as the quantum 
statistical pressure for the Fermion system under consideration. 

e e is the scaling parameter relating the macroscopic coordinates X, T, 

with the microscopic coordinates x,t: X = ex, T = et. The hydrody- 
namic limit is the limit e — ► 0. 

I The embedding from a collection of independent subcubes of periodic 

boundary condition to the cube A E -i. 

I The embedding from a subcube of periodic boundary condition to the 

cube A £ -i . 

u+ £ . . local conservative quantities in a cube of size £ centered at x. 
Ae cube of width £ centered at the origin. 

X s (t,x) A e (i, x) = X(et, ex) 



1.1 A convention 

In the course of our arguments, we will encounter a large (but finite!) number of error terms. 
Therefore, we introduce a common notation that conveys all relevant information about these error 
terms. For k > 1, y any list of symbols, let £l y (x) denote a real-valued function of a; £ M fc , with the 
property that 

lim lim • • • lim £l v (x) = 

where the limits are determined by the names x\, . . . ,Xk of the variables. E.g., any term denoted 
by Q K ,x has the property 

lim limn KjX {£,£) = (1.1) 

£—*oo e^O 

The limits are to be taken in the specified order and the quantities denoted by y = (yi, . . . , yi) are 
kept fixed in the limits. The limit for e is e — > 0, as this is the limit we are considering and similarly 
the limit for £ is unambiguously £ — > oo. The actual value of any Q y (x) may vary from occurrence 
to occurrence. 

For quantities used in error bounds about which no claim of convergence to zero is made, we 
will usually use the notation C y , where y lists the relevant parameters the constant may depend 
on. 
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2 Introduction 



The fundamental laws of non-relativistic microscopic physics are Newton's and Schrddinger equa- 
tions in the classical and the quantum case respectively. These equations are impossible to solve 
for large systems and macroscopic dynamics is therefore modeled by phenomenological equations 
such as the Euler or the Navier-Stokes equations. Although they were derived centuries ago from 
continuum considerations, they are in principle consequences of the microscopic physical laws and 
should be viewed as secondary equations. It was first observed by C. Morrey i n the fifties that 
the Euler equations become 'exact' in the Euler limit, provided that the solutions to the Newton's 
equation are 'locally' in equilibrium. Morrey's original work was far from rigorous and the meaning 
of 'local equilibrium' was not clear. It is nevertheless a very original idea and it contributed signif- 
icantly to the later development of the hydrodynamical limits of interacting particle systems, see 
20 for a review. Instead of considering general classical dynamical systems with two body interac- 
tions, a different approach is to prove as much as possible for some simplified models. Outstanding 
examples are the works by Boldrighini, Dobrushin, and Suhov pP, and Sinai i n the case of one 
space dimension, and the more recent work by Eyink and Spohn who study a <i-dimensional 
classical system of non-interacting particles. In terms of a rigorous proof of Morrey's idea, however, 
significant progress has only been made rather recently ^Hj- This long delay is mostly due to a 
serious lack of tools for analyzing many-body dynamics, in the classical case and even more so in 
the quantum case. 

In this paper, we derive the Euler equations from microscopic quantum dynamics, extending 
the relative entropy method of \2'2\ lltjj to the quantum cases. Our main result was announced 
in [Tl]. As we want to consider the genuine quantum dynamics for a system with short-range 
pair interactions, we cannot take a semiclassical limit. Although one-particle quantum dynamics 
converges to Newtonian dynamics in the semiclassical limit, this limit does not commute with the 
scaling limit needed for the Euler equation. This is most clearly seen in the pressure function, 
for which quantum corrections survive at the macroscopic scale. In fact, one of the conclusions of 
our work is that under rather general conditions, the pressure function is the only place where the 
quantum nature of the underlying system, in particular the particle statistics, survive in the Euler 
limit. 

The Euler equations have traditionally been derived from the Boltzmann equation both in 
the classical case and in the quantum case, see Kadanoff and Baym [S] for the quantum case. 
Since the Boltzmann equation is valid only in very low density regions, these derivations are not 
satisfactory, especially in the quantum case where the relationship between the quantum dynamics 
and the Boltzmann equation is not entirely clear. There were, however, two approaches based 
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directly on quantum dynamics. The first was due to Born and Green who used an early 
version of what was later called the BBGKY hierarchy, together with moment methods and some 
truncation assumptions. A bit later, Irving and Zwanzig jS] used the Wigner equation, moment 
methods and truncations to accomplish a similar result. These two approaches rely essentially on 
the moment method with the Boltzmann equation replaced by the Schrodinger equation. Unlike 
in the Boltzmann case, where one can do asymptotic analysis to justify this approach, it seems 
unlikely that this can be done for the Schrodinger dynamics. 

One of the benefits of our approach is that we develop a general strategy applicable to all 
situations where a number of reasonable assumptions are satisfied. We believe that our general 
assumptions, which are discussed in detail in Section \'2.'dl hold for a large class of physical models. 
We regard proving the properties that we assume as an important, although rather challenging, 
research project in quantum statistical mechanics. 

The main novelty of our work lies in the fact that, for the first time, the relative entropy method 
is applied to a quantum mechanical system. This requires solving a number of technical problems 
which, not surprisingly, all stem from the fact that the local observables corresponding to the 
globally conserved quantities of the dynamics, are represented by non-commuting operators. This 
is mainly discussed in Section |H1 

Although our goal is a derivation of the Euler equations, the relative entropy method indeed 
constructs an approximate solution to the underlying many body dynamics based on solution of 
Euler equations and the concept of local Gibbs states. It thus establishes the key role played by 
the Euler equations: they are not just a set of conservation laws but, with the correct choice for 
the pressure function, they actually dictate the leading approximation to the many body classical 
or quantum dynamics. Thus, the Euler equations may also be used to obtain information about 
the solutions of the many-body Schrodinger equation. 

2.1 Schrodinger and Euler dynamics 

We begin by considering N particles on M 3 , evolving according to the Schrodinger equation 

id t tpt(xi, ■ ■ ■ ,x N ) = Hip t (xi, ■ ■ ■ ,x N ) 
where the Hamiltonian is given by 

N -A- 

^= EY+ E W(xi- Xj ). (2.1) 

j=l l<i<j<N 

Here, W is a two-body short-ranged stable isotropic pair interaction and ipt{xi, ■ ■ ■ , x n) is the wave 
function of particles at time t. We only consider Fermions (such as electrons, but for simplicity we 
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ignore spin) and thus the state space 7i N is the subspace of antisymmetric functions in L 2 (M. 3N ), i.e., 
if)(x ai , — - ,x aN ) = (—l) a ip(xi,--- ,xn), for any permutation a of {1, • ■ ■ , N}. It is more suitable 
not to fix the total number of particles and to use the second quantization terminology. In fact, it 
would be extremely cumbersome to work through all arguments without the second quantization 
formalism. The state space of the particles, called the Fermion Fock space, is the direct sum of 
7i N : H : = 0^ =o H N . 

Define the annihilation and creation operators a x and a x by 

(a x *) N (x u --- ,x N ) = ^NTl^ N+1 (x, Xl ,--- ,x N ) (2.2) 

1 N 

(a+y) N ( Xl ,--- ,x N ) = ^2(-iy-^5(x-x^ N - 1 (x 1 ,... ,£},••• ,x N ), (2.3) 

V 3=1 

where, as usual, ~ means "omit". a x and a x are to be interpreted as operator-valued distribu- 
tions The annihilation operator a x is simply the adjoint of a+ with respect to the standard 
inner product of the Fock space with Lebesgue measure dx. These operators satisfy the canonical 
anticommutation relations 

[a x ,a+] + := a x a+ + a x a+ = 5{x - y), [a+, a+]+ = [a x , a y ] + = , (2.4) 

where 5 is the delta distribution. The derivatives of these distributions with respect to the parameter 
x are denoted by Va x and VaJ. With this notation, we can express the Hamiltonian as H = Hq + V 
where the kinetic energy is given by 

H o = \J Va+Va x dx 

and the potential energy 

V = \ J J dxdyW ( x ~ y) a t a y a v a x ■ 
It is more convenient to put the Schrodinger equation into the operator form, which is sometimes 
called the Schrodinger-Liouville equation. Denote the density matrix of the state at time t by 74. 
Only normal states, which can be represented by density matrices, will be considered in the time 
evolution. Then the Schrodinger equation is equivalent to 

idt'yt = $H7t with 5 H lt ■= [H,jt\- (2.5) 

The conserved quantities of the dynamics are the number of particles , the three components of 
the momentum and the energy. The local densities of these quantities are denoted by u = \i = 
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0, • • • ,4, and are given by the following expressions: 



u x = n x = 



ala 



'X 



1 1 /" 

u x = h x = -Va+Va^ + - / dyW^x - y)a^a^a y a : 



■x 



(2.6) 



We also introduce the notation u = (it 1 , u 2 , u 3 ). This convention will be followed for the rest of the 
paper. We use the bold face for the vector of the conservative quantities and use the frac for the 
vector consisting only the components 1, 2, 3. 

Let Ae denote a cube of width £ centered at the origin. The subscript £ may be omitted if 
it plays no active role. We shall adopt the convention that unbounded observables on A will be 
defined with periodic boundary conditions. E.g., the number of particles in A, the total momentum, 
and the total energy of the particles in A, respectively, are defined by 



In other words, we shall always view A as a three-dimensional torus. 

We slightly generalize the definition of the grand canonical Gibbs states to include a parameter 
for the total momentum of the system: the Lagrange multiplier a. We will work under the assump- 
tion that the temperature and chemical potential are in the one-phase region of the phase diagram 
of the system under consideration such that the thermodynamic limit is unique. The finite volume 
Gibbs states are then given by the following formula: 



The infinite volume Gibbs states (*Jp,a,n are the limiting points of the finite volume ones. It is 
convenient to denote the parameters ((3,a,fi) by A = (A^),^ = 0, • • • ,4 with A = A J = 
(3ai,\ A = (3. Define (notice the sign convention) 




(2.7) 



T r e -/3(H 0iA +VA-aPA.-fiN A ) 



3 




(2.8) 



/u=0 



and 
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These notations allow us to give a compact formula for the unique, translation invariant Gibbs 
state (defined with constant A), as well as for the states describing local equilibrium (defined with 
a>dependent A): 

w A = e |A|(A,u> A/ZA(A) ( 2Q ) 

where Z\(X) is the partition function 

Z A (A) = Tr e l A K A < u >A (2.10) 
The pressure as a function of the constant vector A, is defined by 

V(A) = lim |A|- 1 logZ A (A) 

A^oo 

Denote the expectation value of the conservative quantities in an infinite- volume equilibrium state, 
u>\ introduced following l|2.7j) . by q = (q°, ■ ■ ■ , q ). The we have 

<)r ^aK) (2.11) 



Explicitly, 



1 



p = u) X (n x ) = lim —uj x (N, 



A^R3 |A| 



a; 



e = u; x (h x )= lim 3 ^(H A ) 

Notice that q and e are momentum and energy per volume. 

Again, we will work under the assumption that these parameters stay in the one-phase region, 
the limiting Gibbs state is unique and these definitions are unambiguous. Although momentum is 
preferable as a quantum observable, we also introduce the velocity in order to be able to compare 
with the classical case. The velocity field v (x) has to be defined as a mean velocity of the particles 
in a neighborhood of x. Therefore we have v (x) = q(x)/p(x). We also introduce the energy per 
particle defined by e = e/p. The usual Euler equations are written in terms of p, v, and e. 

In order to derive the Euler equations, we need to perform a rescaling. So we shall put all 
particles in a torus A £ -i of size e~ l and use (X, T) = (ex, et) to denote the macroscopic coordinates. 
For all equations in this paper periodic boundary conditions are implicitly understood. 

The Euler equations for the five conserved quantities, which arise in the limit e — > 0, are given 
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by 




(2.12) 



These equations are in form identical to the classical ones but all physical quantities are computed 
quantum mechanically. In particular, P(e, p) is the thermodynamic pressure computed from quan- 
tum statistical mechanics for the microscopic system. It is a function of X and T only through 
its dependence on e and p. If no velocity dependent forces act between the molecules of the fluid 
under consideration (we consider only a pair potential), the pressure is independent of the velocity. 
The conservative quantities q = (q°, • • • , q 4 ), related to density, momenta and energy as follows: 

q° = p, q l =pv\ q* = e = pe, (2.13) 

In other words q , q 2 , g 3 , and q 4 are momenta and energy per volume instead of per particle as in 
the usual Euler equation 1)2.12(1 . We rewrite the Euler equations in the following form 

^ + E V ^(?)] = ' M = 0,1, 2, 3,4. (2.14) 

4 = 1 

The matrix A is determined by comparison with the Euler equations: 

qi 

5ijP + qiqj/qo (2.15) 

q j (qi + P)/qo ■ 

2.2 Local equilibrium 

To proceed we need a microscopic description of local equilibrium and a microscopic prescription 
to compute the pressure from quantum statistical mechanics. Suppose we are given macroscopic 
functions q(-X"). We wish to find a local Gibbs state with the conserved quantities given by q(X). 
The local Gibbs states are states locally in equilibrium. In other words, in a microscopic neigh- 
borhood of any point x S T 3 the state is given by a Gibbs state. More precisely, we wish to find 
a local Gibbs state with the expected values of the energy, momentum, and particle number per 
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unit volume at X given by q(X). To achieve this, we only have to adjust the parameter A at every 
point X. More precisely, we choose X(X) such that the equation (|2.1H) holds at every point, i.e., 

dx^(x) - q {x) - 

If we denote the solution to the Euler equation by q(X,T), then we can choose in a similar way a 
local Gibbs state with given conserved quantities at the time T. Define the local Gibbs state 

ul = -^exp [e- 3 (A(et,e-), u^.J (2.16) 

where c e (t) is the normalization constant. Clearly, we have that uf{ux) = q^(ex,et) to leading 
order in e. Our construction of local equilibrium states is consistent with the abstract framework 
discussed in [2"T] . 

Later we will need the following relation for the normalization constant c e (t): 

^logc e (t) =eTr uf[e- 3 {d t \(et, £•), u}^] (2.17) 

Since the inner product almost exclusively taken on A E -i, we shall drop this subscription or replaced 
it by ( , ) e -i for the rest of this paper. 

In summary, the goal is to show that, in the limit e — ► 0, the following diagram commutes: 

q(X ,0) q(X,T) 



local equilibrium 



limit £ J. of expecta- 
tion of locally averaged 
observables 



Schrodinger 

70 * 7e-lT 

As smooth solutions of the Euler equations are guaranteed to exist only up to a finite time ^U] , 
say To, we will formulate our assumptions on the dynamics of the microscopic system for a finite 
time interval as well, say t £ [0, Tq/e]. Note the cutoff assumptions below would hold automatically 
for lattice models. 

2.3 Assumptions and the main theorem 

Our main result is stated in Theorem 12.11 below. First, we state the assumptions of the theorem 
with some brief comments. There are three kinds of assumptions. 

The first category of assumptions could be called physical assumptions on the solution of the 
Euler equations that we would like to obtain as a scaling limit of the underlying dynamics, and on 
the pair interaction potential of this system. 
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I. One-phase regime: We assume that the pair potential, W, is C 1 radial and supported in a 
ball of radius R. Furthermore, we assume that W is stable in the sense that 

W(x) = W (x) + Wi{x) where W > , W (0) > and W x is positive definite. (2.18) 

Here, the positive-definiteness of W\ refers too the sesquilinear form, not the function itself. I.e., 

W\{xi — Xj)z~iZj, for all, n > 1, Xj, Xj G R 3 , Zi, Zj G C. 

l<i,j<n 

Such potentials automatically satisfy the usual super-stability property ^Sj. In particular, they are 
stable, i.e., there is a constant B > such that, for all N > 2, x\, . . . , xn G K 3 , 

W{ Xi - Xj ) > -BN 

i<i<i<^ 

Of the Fermion system with potential W we assume that there is an open region D C M 2 , 
which we will call the one-phase region, such that the system has a unique limiting Gibbs state 
and a regular pressure function for all values of particle density and energy density (p, e) G D. 

The solution of the Euler equations we consider, q(X,T), will be assumed to C 1 in X for T £ 
[0,To], and have local particle and energy density in the one-phase region for all times T € [0,To]. 

I. e., {p(X,T),e(X,T}) G D, for all X G Ai and T G [0,T ]. 

The next category of assumptions is on the local equilibrium states for the Fermion system that 
we construct and on their time-evolution under the Schrodinger equation. 

II. Cutoff assumptions: Suppose that jt is the solution to the Schrodinger equation ()2.5|) with a 
local equilibrium state as initial condition, constructed with the parameters derived from a solution 
of the Euler equations (with the appropriate pressure function) for times t in a finite interval [0, To], 
that does not leave the one-phase region. We make the following two assumptions. 

1. High-momentum cutoff assumption: Let N p (t) = Tr 7ta+ a p , where ajf is the Fourier transform 
of ax - Then there is a constant c > such that for all t < Tq/e, 

e d j dpe cp2 N p (t) < C To (2.19) 

where Ct , is constant only depending on To. 

2. Non-implosion assumption: There is a constant Ct (not necessarily the same as Ct in the 
previous paragraph) such that for all t < Tq/e 

1 2 1 

< C To (2.20) 











Tr < 


j t s d / dxn x 


/ n y dy 








J\x-y\<2R 
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where R is the range of the interaction W . 

Finally, we have an assumption on the set of the time-invariant ergodic states of the Fermion 
system. To state this assumption we need the notion of relative entropy, of a normal state 7 with 
respect to another normal state u. Let 7 and u> denote the density matrices of these states. The 
relative entropy, £(7 | u), is defined by 



For a pair of translation invariant locally normal states, one can show existence of the relative 
entropy density [T^], defined by the limit 

s(j\u) = lime 3 5(7A | omJ , 

where 7a _ x and lo\ _ t denote the density matrices of the normal states obtained by restricting 7 
and u> to the observables localized in A e -i = e~ 1 Ai. The existence of the limit can be proved under 
more general conditions on the finite volumes, but this is unimportant for us. 

III. Ergodicity assumption ("Boltzmann Hypothesis"): All translation invariant, ergodic 
with respect to space translations, stationary (i.e., time invariant) states to the Schrodinger equation 
with the Hamiltonian H are Gibbs states with the same Hamiltonian provided they satisfy the 
following assumptions: 1) the density and energy is in one phase region. 2) The relative entropy 
density with respect to some Gibbs state is finite. 

We expect that the cutoff assumptions hold for the solutions 7t of the Schrodinger equation 
that we employ, but for now there is no complete proof that it holds for Gibbs states other than 
the free Fermi gas. For Gibbs states in the high temperature region we expect these assumptions 
can be proved by using some type of cluster expansion methods. A partial result in this direction 
has been obtained recently, in the case of Bosons, by Gallavotti, Lebowitz, and Mastropietro in [B]. 

For the rest of this paper, we shall assume this cutoff assumptions for the solution to the 
Schrodinger equations as well as the Gibbs states in the one phase regions considered in this paper. 

We wish to point out that in the treatment of the classical case in jlbj the cut-off assumption 
2 was not needed. There is however no proof for the cut-off assumption 1 even in the classical 
case. (In the usual quadratic kinetic energy was replaced by one with bounded derivatives 
with respect to momentum. So the cut-off assumption 1 is not needed too.) 

The cutoff assumptions are technical in nature. For Fermion models on a lattice instead of 
in the continuum, no cut-off assumptions are required. The Boltzmann hypothesis on the other 
hand is a fundamental problem in statistical physics. A version of it was proved to hold for a 




Tr {7(log7 — logw)} if kerw C ker7 



otherwise 
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classical ideal gas by Eyink and Spohn in [5]. Gurevich and Suhov f7] proved that a stationary 
Gibbs state to a classical dynamics with a Hamiltonian H has to be a Gibbs state with the same 
Hamiltonian. Under the assumption that the stationary measures velocity distribution has no 
correlation (a weaker assumption than in [7]), the Boltzmann hypothesis was proved for classical 
gas with two-body interaction [16,. 

Our main result is the following Theorem. We also expect it to hold for Bosons with a super- 
stable interaction. 

Theorem 2.1 Suppose that q(X, T) is a smooth solution to the Euler equation in one phase region 
up to time T < Tq. Let lo\ he the local Gibbs state with conserved quantities given by q(X, T). 
Suppose that the cutoff assumptions and the ergodicity assumption hold. Let "ft be the solution to 
the Schrddinger equation (j2,5() and 70 = ujq (Note that jt depends on e). Then we have 

lim sup s(7t|wf) = 

£ -*°0<t<e- 1 T 

Ln other words, uf is a solution to the Schrddinger equation (|2.5|) in entropy sense. In particular, 
for any smooth function f on A, we have, for all < T < Tq, 



lime 3 / dxf{ex) [j e -i T {vL x ) - q(T,ex)\ = 

*-*0 A _i 



The proof of this theorem will be given in Section 17.21 
Illustrated with a diagram, the main theorem says 



Euler equation 

<0 > 



lim e i0s(7 e -i T l^q T )=0 

Schrddinger equation 

70 > 7e-iT 

Notice that we have proved more than just convergence to the Euler equation. We have shown that 
the local-equilibrium Gibbs state constructed from the evolution of the Euler equations solves the 
many-body Schrddinger equation, approximately in entropy sense. 

2.4 Outline of the proof 

The basic structure of our proof follows the relative entropy approach of [161 122j . The aim is to 
derive a differential inequality for the relative entropy between the solution to the Schrddinger 
equation and a time-dependent local Gibbs state constructed to reproduce the solution of the Euler 
equations. The time derivative of the relative entropy can be expressed as an expectation of the 
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local currents with respect to the solution to the Schrodinger equation. Since we do not know the 
solution well-enough, this expectation can not be computed. 

Step 1: Replace the local microscopic currents by macroscopic currents. The basic idea in 
hydrodynamical limit is first to show that the local space time average of the solution is time 
invariant. From the Boltzmann hypothesis, ergodic time invariant states are Gibbs. For Gibbs 
states, we can replace the local microscopic currents by macroscopic currents. This is the first step. 
In the quantum setting, there are several crucial issues we need to address. 

la: Construct a commuting version of the local conserved quantities. Recall macroscopic cur- 
rents are functions of the local conserved quantities, i.e., density, momentum and energy. For 
the microscopic quantum system, the local conservative quantities are operators which commute 
only up to boundary terms. In order to express the macroscopic currents as functions of the local 
conserved quantities, we need either to prove that the non-commutativity does not affect the macro- 
scopic currents or we need to construct some commuting version of the local conserved quantities. 
As the first approach seems very difficult to carry out, we follow the second one and construct a 
commuting version of local conservative quantities in section 4. 

lb: Restriction to the one phase region. Since the Boltzmann hypothesis holds only in the 
one phase region, we have to exclude the region outside the one phase region. To perform this 
restriction to the one-phase region , we would normally multiply the observables by some cutoff 
function. In our case however, the cutoff function does not commute with the local currents. This 
seemingly trivial multiplication by a cutoff function illustrates the kind of technical problems we 
have to address in this work. Our approach to this is presented in Section El 

lc: Virial Theorem. Even assuming the local ergodic states are Gibbs in the one phase region, 
we still have to compute the macroscopic currents from the microscopic currents. This requires a 
virial Theorem, which we provide in Section ITU1 

Step 2: Estimate all errors by local conservative quantities. As will become clear, the errors 
associated with the cutoff of the one phase region are difficult to control directly. We shall bound 
them by the local conserved quantities. This will be carried out in sections 5 and 6. 

Step 3: Derive a differential inequality of the entropy with error term given by a large deviation 
formula. After Step 2, we have an expression of the derivative of the entropy in terms of local 
(commuting) conservative quantities. Since these quantities commute, by an entropy inequality, we 
can bound it by a large deviation expression. Notice that it is crucial that we control everything 
by commuting objects. There is no large deviation theory for non-commuting observables. 

After this step, the standard relative entropy method provides the rest of the argument. Tech- 
nically speaking, the main difficulty to study a quantum mechanical system, in comparison with 
a classical one, can be traced back to the non-commutativity of the algebra of observables. E.g., 
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suppose A and B are two self-adjoint operators representing observables of the system. A simple 
inequality, such as |vl + -B| < |A| + |£>|, which is used numerous times in estimates for classical sys- 
tems, is false, and so is \AB\ < \A\ \B\. Therefore, there are essentially no absolute values taken in 
our proof and we estimate all quantities by commuting versions of the locally conserved quantities. 
Of course, these inequalities hold with the absolute value replaced by the norm. However, we will 
frequently deal with error terms that are expectations of unbounded observables, such as, e.g., the 
high-momentum contributions to the energy. Clearly, norm estimates are useless in this situation. 



3 Relative entropy identity and high momentum cutoff 
3.1 Entropy identity 

The first step in the derivation of a diffential inequality for the relative entropy is the compute the 
derivative. Suppose 74 is a solution to the Schrodinger equation. Recall that one has 

-^Tr A(t)B(t) = Tr A'(t)B(t) + A(t)B'(t), -^Tr e A W = Tr e A ®A'{t) (3.1) 



and 



d 
dt 



5(7*) = 



Thus we have for any time-dependent density matrix pt the identity 



(3.2) 



4; 5 (7t I = Tr 1/t{-i$H ~ d t }log p t 
dt 



(3.3) 



This identity replaces the relative entropy inequality in [23|25|. Thus, by (|2.16j) . 



d_ 

dt 



s{ lt \uj £ t ) =Tr 7t {-i5 H - d t } {(\ £ {t,-),u) £ -! - e 3 log c £ (t)} 



(3.4) 



where ( , ) £ -i = ( , )a £ _ 1 . 

Let w denote the current tensor with components w£ x defined by 



w k,x 



W 



w lx 



■2'* 



W\x - y) 



(x-y)j(x-y) k 



\x - y\ 



CL X CLy CLyCL x , 



k = 1,2,3 



1 

A L 

i 



Vfca+Aa^ - Aa^V k a x 

(x - y) k (x - y)j 



(3.5) 



(3.6) 



+ - I dyW(x - y) V k a x Oy a y a x - a x ciy a y V k a x (3.7) 



dy 



W'{x-y)- 



f - y\ 



a x Vja^a y a x - a x ayVja y a x 



(3.8) 
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where we have used the rotation invariance of the potential to write 

„,„ x dW(r) l 

We have the following proposition. See Section for the derivation of these expressions for the 
current. 

Proposition 3.1 Let Q,\{e) be defined by the equation 

3 

i5 H (X £ (t, -),u) £ -i =eJ^(V i A e (t,-), Wj(t)) e -i + n x (s) (3.9) 

Then Q\(e) is an error term which satisfies the condition 

limTr 7^ A (e) = 

£^0 

These expressions of the microscopic currents seemingly bear no relationship to the macroscopic 
currents in the Euler equations, even when one assumes that 7$ is locally Gibbs. This difficulty 
already appears in the classical case. But by reasonably straightforward computation and appli- 
cation of a quantum version of the virial theorem 111). II one can show that indeed these currents 
correspond to the standard Euler equations given in (|2.12|) . 

Define Vo = dt and Wq x = u%. We have 

3 

{-iS„-d t }e 3 (X £ {t,-),u) £ -i =-£^(V,-A £ (V), Wj -(t)) E -i :=eG(X £ ,a + ,a) (3.10) 

3=0 

where \ £ (t,x) = \(t,x). Introduce the notations V = (Vo, V), w = (wo, w) and 

3=0 M=0 j=0 

Then we can rewrite the last expression as 

{-i6 H -d t }e 3 {X £ (t,-),u) E -i =-(VA £ (t,-) • w(t)) E -i (3.12) 

If we wish to emphasize the dependence on the operator, we shall write Wj x = Wj x (a + ,a). From 
now on, we shall drop the subscript e _1 in ( , ) e -i. 
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3.2 High-momentum cutoff 

Most of the estimates we need are obtained using bounded versions of the creation and annihilation 
operators, i.e., suitable so-called smeared operators. Physically, this corresponds to introducing a 
high-momentum cutoff. The precise form of the cutoff will be important for us, as we will have 
strict requirements on the behavior of the error terms for the proof to go through. 
Let (pM be a smooth function such that 

1. \<f> M (p) ~ 1| < e~ M2 for \p\ < M and \4> M (p)\ < e~ M2 for \p\ > 2M. 

2. The support of 4>m is bounded in a ball of radius e M2 . 

To construct such a function, let g be a smooth function supported in |x| < 2 such that 

g(p)<C[l+ P 2 ]- 3 
Define g\(x) = g(x/X)X~ 3 ^ 2 . Let hu be a smooth function such that 

h M (p) = 1 for \p\ < M and h M (p) = 0, for \p\ > 2M . 

Let 

4>m = (g\ *g\)hM 

Notice that f 4>m = 1- Let A = e M2 . Then we can check easily the properties 1 and 2. Although 
4>M is supported in a ball of radius e M2 , its mass is concentrated in a ball of radius M -1 . More 
precisely, there is a constant c such that 



/ h,M{x)dx < e 

J\x\>r 



-crM 



Define 



= J 4>m{x- y)a+ = a+ m , a xM = j <t>{x - y)a y = 



<t>x,M 



where 4> x ,m = 4>m{x - y). In our setting 4>x,M = 4>x,M- Define 

Notice that Va^ M and M are bounded operators localized in a ball of radius e M2 . 

We now perform the preliminary truncation. Denote the cutoff version of the current by 

v?j,x,M = w jjX (alj, a M ) (3.13) 

Notice that vfj, x ,M is bounded. The difference between the kinetic energy term in the energy current 
with and without cut-off can be calculated using 

Va+Aa x - Va+ M Aa Xi M = Vb+ M Aa x + Va+Ab^M 
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where 



bt 



a' — a 



xM 



Lemma 3.2 For any state 7 that satisfies the cutoff assumptions and for G defined by (|3.1Uj) . we 

have 



e d Ti 7 / dx 



Vb+ M Aa x + VatAb x 



'xM 



< e 



-cM 2 



Proof: By using the Fourier transform, we have 

J dxVb+ M Aa x = J dp(l - cf>(p))p 2 a+ a p 
Let N p (t) = Tr 7a+a p . Then 



Tr 7 J dp(l - <fi M (p))p 2 a+a p = J dp(l - (j>M{p))p 2 N p 

The lemma is thus a simple consequence of the Chebeshev inequality and the definition of 4>M- 
From the Schwarz inequality 

Tr je d J dxVa^ J dyW {x — y)a+ a y a x 



1/2 r 

< e a <! Tr 7 / dxVa^Va x } I Tr 7a 



dyW(x - y)n y 



1/2 



< C 



Thus for any state 7 satisfy the cutoff assumptions, we have 



Va+ay-a y a x - Va+ M a+ u a yM a x 



x,M yM V' 1 



Tr je d J dxdyW(x - y) 
Lemma 3.3 For any state 7 satisfy the cutoff assumptions, we have 

Tr 7 e o! [G(A,a + ,a) -G{X,dt ,a M )} < e~ cM2 



M 



< e 



-cM 2 



Thus we have 



^s(7iK £ ) = ^Tr j t G{\ E ,a^,a M ) - £ 3 ^ log c £ (t) + £ l M 



(3.14) 



(3.15) 



(3.16) 



with £\{ < Ce~ cM2 . The precise form of the last estimate is crucial as we shall see later on. 

We recall the crucial relative entropy inequality |15j : for all self-adjoint observables h, and for 
any 5 > 0, : 

j(h) ^(TMogTr e 5h+losuJ + 8~ 1 S(j\lo) (3.17) 
A proof will be given in the appendix. 



21 



4 Construction of local commuting observables 



The conserved quantities commute as global observables on A with periodic boundary conditions. 
In fact, this is essential for the classical equations of motion to make sense. For any bounded 
quasi-local observable X on the Fock space define 

X(q) = Tr u x X (4.1) 

where the chemical potential A is the dual of q in the sense of (|2.11|) . We define X only for 
arguments in the one phase region. 

Local averages of the densities of the conserved quantities, however, do not commute due to 
boundary effects. Therefore, we cannot extend the functions X to functions of the operator-valued 
local densities of the conserved quantities, which is what we would like to do. To circumvent this 
difficulty, in this section we construct commuting versions of the local conserved quantities. 

We have for any smooth function J 

Tr 7t( J(et, e ■ ■ ■ ), X) = Tr 7t Av J(et, ex) Av t z X + QAe) 

zeA £ -i \z-x\<£/2 

where Av(-), stands for the average of its argument over the domain indicated in the subscript, and 
lim e ^o = f° r an Y fixed. 

Denote := ua^ := t x udx the local conservative quantities and we would like to replace 
the microscopic current Av t z X by certain function of the local conservative quantities t z U£. 

\z-x\<£/2 

Unfortunately the components of do not commute and functions of are not well-defined. 
In fact, even the definition of is ambiguous since we did not specify the boundary condition. 
Intuitively, the components of actually commute up to boundary terms, and the ambiguity should 
be negligible in the limit i — > 0. Since it is rather difficult to control these boundary terms in a 
simple way, we construct in the following a commuting version of the local conservative quantities. 

4.1 Construction of an isometric embedding 

Let / be a smooth function with 

f(s) = 1/V2 ifs<0 =0 if s > 1 

and 

/(I) = /'(!) = /"(I) =0. 

For any given rj, < rj < 1/2, let 

1/2 

g(t) = 1 - f I ^— I , < t < 1/2 
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g(t) = f (^i-j , t > 1/2, 

g(t)=g(-t), t€R. 
Then g is smooth, supported in \t\ < 1/2 + rj and 

Let x(x) = g(x 1 )g(x 2 )g(x 3 ). Then 

j> 2 (*+j) = i ( 4 - 2 ) 

Let := T a Af be a cube of size £ ± 4£r? centered at a. Let 

Xa(af) = X((x - a)/t) 

be a smooth function supported in r Q A^ + 2&j C a + and Xa( x ) = 1 in r a A^_2^ D a - . We collect 
these relations in the following: 

a~ = TaAg-Mr, C T a A t _ Ur) C {x : Xa(«) = 1} 
C {x : Xa(x) / 0} C T a A i+2 e v C T a A t+Urj = a + (4.3) 

There is a wide range of choices for rj. The main restrictions we needed are 

r]£ — > 00, 77 — > . 

We shall choose, for simplicity of notation 

v = r 1 / 2 

for the rest of this paper. 

Recall the configuration space 5(A) is the space 

5(A) = {x := (xi, ■ ■ ■ , x n ) : n £ {0} U N, xj £ A for all j} 

Denote by T(A) the space of antisymmetric functions from the configuration space 5(A) to the 
complex numbers. With the standard L 2 inner product, T(A) is a Hilbert space. 
Define I A from T(a+) to T(A) by 



(7 A ^)(xi,--- ,x n ) 



YlXa(Xj) 



ip(x 1 , ■■■ ,X n ) 
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Usually we shall take A = A e -i and omit the labels A and a whenever they are obvious or unim- 



portant. It is crucial that 



is symmetric w.r.t. permutations of x so that I^ip is 
antisymmetric as a function of x. Define /* to be the adjoint of /, i.e., we have (/*/, g) = (/, Ig). 
Let X be an observable X on a defined by 

X= dx 1 ---dx k dy 1 ---dy k f(x 1 ,...,x k ;y 1 ,...,y k )a^ ua ---a+ k;a ay kta ---a yua . 

where / is a distribution (kernel) with support in (a + ) x2fc . Here we have labelled the operators by 
a to emphasized the cube a. We can check the identity: 



I* XI = J dx 1 ---dx k dy 1 ---dy k Xa(x 1 )---Xa(x k )xa(yi)---Xa(yk) 



xf(x 1 ,...,x k ;y 1 ,...,y k )aZ 1 ---aZ,a 



(4.4) 



*x fc "'J/fc "W 

as an operator on the torus A e -i. 

From this definition the pull-backs of all observables we need, i.e., I* XI for a conserved quantity 
or current given by X, can easily be computed. By using the appropriate distribution kernels /, 
observables involving derivatives are included. E.g., 



I* 



/ dxf(x)Va+a x 

Ja+ 



For the kinetic energy we have 



f dxf(x)Vji 



I* - dy5'(y - x) dxf{x)a+a x 
-J dy5'(y-x) J dxf '(x)x a (x)x a (y)ay a x (4.5) 
J dxf(x) [xa(x) 2 Va+a x + Xa(x)Vxa(x)a+a x ] 



dxf(x) [x a {x) 2 Vja+Vja x + (V ' jXaix)) 2 a+ a x (4.6) 
+Xa(x)V j x a (x)(V j a^a x + a+Vj-a^)] 



If we take / = l a -(x), we have f(x)x a (x) = 0. Together with Xa(x) = 1 in r a A^_ 2 ^ D a , we 
have 



I* 



L 



dxVaZa x 



dx\7ja x Vja x 



L 



dx [Va^a a 



dxV ja^V jd x 



(4.7) 



4.2 Commuting local conserved quantities 



Let H a + , P a + be the total energy and momentum operators on a + with periodic boundary condition. 
Then H a + and P a + commute with each other and also with the number operator N a + . Denote 

u Q+ =r 3 (P a+ ,N a+ ,H a+ ) 
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Since the image of I is in the domain of H a + and P a +, the operator I*u a +I is well-defined. Since 
the components of u a + commute, the function 

I*X(u a+ )I 

is now well-defined. We shall use the notation 

u te : = u «+' n ti : = ^ N a+, h+ £ := r 3 H a+ (4.8) 
when a is centered at x. When x = 0, we shall omit the subscript x. 

4.3 Local average of currents 

Let 

w„± = / dxw x , W a ± = / dxw x 

be the average over the cube of the currents w x , where we have divided the integration by £ 3 
which is approximately the volume to the cube . by definition mv x is an operator on the torus 
A e -i. Since a| can be viewed as an operator on a + with periodic boundary condition for x G a + , 
w Q ± can be understood as an operator on the cube a + as well. We shall use the same symbol in 
both contexts. 

Recall the cutoff version of the current ~w x ,M = w x (a^,aAf) (j3.13() . Thus we can define the 
cutoff version of the current 

w m,q± = / dxw M ,x (4.9) 

Here ~w M ^ a ± can be viewed as an operator either on the torus A e -i or a + with periodic boundary 
condition. 

Recall that ajj, M are bounded operators localized in a ball of radius e M2 centered at x. Thus 
the support of ajj. M is contained in r a A^_ 2v ^ for x £ a~ = T a A e _ 4 ^j, as long as \/l > e M2 , which 
we shall assume from now on. Since Xa{x) = 1 for x £ T a A £ __ 2 ^, following the proof of (|4.7j) we 
have the identity 

J*w M)Q -I = w MjQ -; (4.10) 

here w Ma - is understood as an operator on the torus A e -i on the right side and as an operator 
on a + on the left side. Define the notation 

where K,t = T * A e-4Vl = T * A i- 
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From l|4.10|) . the boundedness of a M x and simple counting of the number of terms, we have the 
following lemma. 



Lemma 4.1 For any state 7, and any smooth function J, we have 



Tr 7 <J(e.),w M > = Tr 7 AvJ £ (x)I*w 



M,x,£ 



j + n M (e,£) 



(4.11) 



where the error term Qm(£,£) vanishes in the sense given by (jl.l|) and 

Av = e 3 / dx 




Applying this Lemma to a smooth function J(et, ex) and average over t, we have 



Av Tr -y t (J(e;et),w M ) = Av 7t ( J e (t, ■) , I*w^ , t I) + n M (e, £). 



(4.12) 



5 Bounds on the currents 

The aim of this section is to show that the currents with momentum cutoff, i.e., the quantities 
^A/o±) can De bounded by a multiple of the Hamiltonian plus particle number. This is the 
content of the following lemma. 

Lemma 5.1 The following operator inequalities hold: 



on the left side is replaced by W m a - . 

As it is essential for the proof of this lemma, we first recall a standard stability result based on 
the superstability conditions we have assumed on the potential W. Stated in words, the result says 
that a large class of two-body quantities can be bounded in terms of the two-body interaction and 
the particle number. We state this result as a lemma for functions but it obviously extends to the 
corresponding second quantized observables. 

Lemma 5.2 Suppose U is a positive bounded function with compact support on IR 3 and W is a 
superstable potential stated in the sense of (j2.18|) . Then there is a 5 > such that 



W M , a+ <CM[H a+ +N a+ ] 




Note that the dependence on M in the right hand side is linear. Similar inequality holds if a + 



N N 



5 U(x a -x p )<Y, W(x a - xp ) + N . 
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The proof of Lemma 15.21 is contained in Ruelle's book ^Sj or see |16j . 
Proof of Lemma I5.lt First, we treat the one-particle (i.e., quadratic in the a^'s) terms and 
show that they can be bounded by the kinetic energy term and chemical potential term of the 
Hamiltonian. 

We will use the following inequalities several times without further reference: for any pair of 
bounded operators A and B, and c > 0, one has 

A*B + B*A < cA*A + c~ 1 B*B 

and 

A + B<\A\ + \B\ . 

From the last inequality it follows that for a bounded family of self-adjoint operators A x , and a 
real valued L 1 -function /, one has 



f(x)A x dx< J \f(x)\\A a 



\dx 



Note that one cannot replace the LHS of these inequalities by their absolute values unless all terms 
commute. 

We start with bounding the momentum components of ~Wm : 



< 



/ dx V k a+ M a XjM - a+ M V fc a Xj M 
Ja+ 1 ' J 

/ dxV k a+ M V k a x ,M + 0,+ M a X! M (5.2) 



We would like to obtain bounds by multiples the Hamiltonian and the number operator without 
momentum cut-off M. For this we use the following inequalities: 

a t,M a x,M < C\(f> M \*a x a x (5.3) 

Vo+ M Va X)Af < CM\V((>M\*a+a x (5.4) 

Va x,M Va z,M < I^Afl * Va+Voa, (5.5) 

Aa+ M Aa x , M < CM\V4> M \ * Va+Va x (5.6) 

In the above expressions the convolutations are with respect to the variable x in the RHS. The four 
inequalities are proved in almost identical fashion. E.g., the first inequality is obtained as follows: 



a x,M ax > M = J dzdw(f) M {x - z)aj4> M (x - w)a w 

< - j dz [ dw\cf>M (x - z)\\((>m (x - w)\ [a+a z + a+a„ 



2 . 

< \(Pm\ * a x a x 
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where we have used that J \(f>M\ = 1- For (|5.4|) and ()5.6|) one also has to use ||V0m||i < CM, for 
a suitable constant C, but otherwise the proofs are the same. Now, we can finish the bound of 
Mai by using ((OJ) and in (1Tk2J). We obtain 



^fc,M,a+ ^ J + dx\(j) M \{x-y) [V k a+V k a y + 

< C(tf 0) a + iV Q ) 



By stability of the potential the kinetic energy term in this bound can be replaced by the full 
Hamiltonian, up to an adjustment to the constant C. This completes the proof of the lemma for 

For the other one-particle terms of Wj\/ jQ one proceeds in the same way. E.g., for the last 
inequality one starts from 



< MV k a+ M V k a xM + M _1 Aa+ 



The rest of the argument is the same. In summary, the results are 



dx 



VfcCi M a xM — 0- M^kO> x ,M 



dx 



< C(H a+ +N a+ ) 

/ dxVjdt M V k a XiM < C(H a ++N a +) 

Ja+ 

V k a+ M Aa x , M ~ Aa+ M V k a x , M < CM(H a + + N a - 



The two-particle terms (quartic in the aft's) appearing in ()3.6j) and ()3.8|) . we can follow the 
same procedure. E.g., to bound the middle term of 1)3. 8|) . we start from 



^ k a t } M a y,M a y, Ma x,R4 ~ a t,M a y,M a yM^ k&x,M 



(5.7) 



- Ma t,M a y,M a vM a 3:,M + M ^ a t,M a y,M a V,M^ a x,M 
The first term can further be bounded in a way similar to 1)5.3)1 : 

a x~,M a y,M a y> Max ' M - J dz \<pM(z - y)\a+ M aya y a x ^ M 
= J du(j)M(u - y)a+a+ M a XtM a y 
< J dudv\(f) m \(u - y)\(f) m \(v - x)a x ~a y v a y a x . 
The same quantity appears in 1)3.6 J) . In both cases, after integration, we get something of the form: 
Ml dx j dy[ dudv\cf)M\(u — x)G(u — v)\(/)m\(v — y)]a x a,ya y a x , 
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where G is a non-negative function of compact support. To estimate this term, we use Lemma 15.21 
to obtain a bound of the last expression of the form CM times the potential energy. 
The second term in the RHS of (|5.7|) gives rise to 

M _1 / dx j dyM[ dudv\V(t>M\{u — x)G(u — v)\(/)m\{v — y)\a£ at a y a x 



and something of the same form for the third term in which, with another application of 

Lemma l5.21 can also be bounded by the CM times the potential energy. This completes the proof 
of the lemma. ■ 



6 Local ergodicity 

Recall that by assumption the solution up to time t < Tq/e of the Euler equations has density 
and energy taking values in a compact set strictly contained in the one phase region of the phase 
diagram of the fermion systems. Let a K be a smooth function supported in the one phase region 
such that ct k = 1 on this compact set. Furthermore, we require that as k — > 0, cr K becomes the 
characteristic function of a compact neighborhood of this set contained in the one phase region. 
Since the phase transition region depends only on the density and energy, a K needs to depend only 
on the density and energy. We will take <r K (e, n) of the form erf (e)o"2 (p) wnere of and are some 
smoothed versions of the characteristic functions on a set of sufficiently high e and sufficiently low 
p, respectively. 

The aim of this section is to prove the following theorem. Recall X is defined in Q4.1|) . 
Theorem 6.1 For all smooth functions J, and X any one of the components of w A/ we have 
Av y t AvJ(et,£x)I*X-J 

t<T/e % 

< Av Tr 7( AvJ(rt, £ .) {f(W)(u+ f )7 + I* (1 - ^(u+^pQ^l - a"(u+))J } 
+n.j^x(e,£,a) (6.1) 

where a K = y/a K (2 — a K ). 

The function a K behaves essentially the same way as a K , i.e., it is a smooth version of a char- 
acteristic function supported in the one-phase region. 

As a first step towards the proof of Theorem 16.11 we partition A e -i into cubes of size ae _1 , 
where a is a sufficiently small positive constant. For any z S A e -i, let Q = A z ae ~i denote the cube 
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of size as 1 centered at z. For any bounded quasi-local observable Z, define the average of Z in 
the cube Q by 



Z Q = Avr^Z 



We also divide the time interval [0,e _1 T], into disjoint intervals of size 2ae~ 1 and label the centers 
by tx, ■ ■ ■ t n , n = T I (2a) (the n-th interval is [t n — a, t n + a] n [0, e _1 T]). 
Since J is a smooth function, 



(J e (t, •), Z) = rC 1 £ Av [ J(etj,ez) { Av Tr lt Z K \ 



+ ^ Zi j(a,e) (6.2) 



where lim a _>o lim £ _>o ^Z,j{ a j e) = and the average is over z £ ae 1 Z 3 n A e -i. 

For Q, a,j fixed, define a family of states labelled by e consisting of the states defined by 

tF(^) = , Av Tr lt Z Q 
\t-tj\<a/e 

Then | e > 0} is w*-precompact and, hence, has at least one limit point. 

Lemma 6.2 Let uj be the Gibbs state on A e -i defined in (|2.9f) with A = A e -i and the chemical 
potential A := A is chosen to be 

A = Av x X(0, ex) 

where A(0, •) are the parameters for the initial condition defined in l|2.16|) . Then for any t > the 
relative entropy 

sint I w) <c 

/or some constant C depending only on the initial value A(0, •). 
Proof: Recall the initial state is 

w o = ^Lexp[ £ - 3 (A £ (0,0,u)] (6.3) 

Then we have 

s(wg|w) = y ttewg f (A £ (0, •) , u) - (A , u)) + e 3 log c e (0) - e 3 log Za,^ (6.4) 

Since each term on the right side is bounded, we have s(cjq|u)) < C . From a simple direct calculation, 
we know that s{^\uj) is a constant of motion. This proves the Lemma. | 
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Lemma 6.3 Fix the parameter a and let n be any limit point of {7?'"' | e > 0}. Then n is 
a translation invariant, time invariant state of the dynamics. Furthermore, the specific relative 
entropy of r\ with respect to the translation invariant state uj-^, satisfies the bound 

Proof: The invariance under space and time translations is an immediate consequence of the 
scaling by e . Since the proof for quantum case is parallel to that of the classical case, we refer 
the reader to JS] for a proof of the classical case. To show that the specific relative entropy with 
respect to uj^ is finite, we start form Lemma 16.21 stating that the relative entropy e 3 S(^ \ uj^) < C 
for a suitable constant C. 

The operations of averaging over translations in a cube Q and over times in an interval [t{ — 
a/e,ti +a/e], are completely positive, therefore, by the monotonicity (or convexity) of the relative 
entropy (see, e.g., US]), we have 

e 3 S( Av Av 7 f ot v \u£)<C 

\t-tj<a/eyeQ 

The relative entropy is also monotone with respect to restriction to the algebra of observables of a 
subvolume. Therefore we have 



e Si Av Avt^otJ 

It-tjKa/ey&Q y ' 



<4Iq) ^ c 



Now, 77 is a limiting point of {7^ ' 3 \ e > 0}, where 

ryQ,j = Av Av 7? o Ty 
\t-tj\<a/eyeQ 

Therefore, by the lower semicontinuity of the specific relative entropy, we can conclude 

< (2a)- 3 limsupe 3 5(7 e Qj ' | w|) < C(2a)- 3 

e 

■ 

Consider any limiting point n of {^' J \ s > 0}. Since n is translation invariant, we can 
decompose it into ergodic components (with respect to space translations) and there is a probability 
measure fi supported on ergodic states u such that 



n = j u) fi(dui) 
The key property of rj is the following lemma. 
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Lemma 6.4 Let n be as above, and X £ ^4a . Then there is £1m(£, k ) such that 



r,(I*X7I)-r, 



f(W)(u+)/ + r(l-a K (u+))X £ -(l-a K (u+))/} 



<n M (i,K) (6.5) 



where a K = y / cj k (2 — <t k ). 

Proof of Theorem 16.11 assuming Lemma 16.41 

Let 



Z = I*X7I 



r(d K Xa K )(u+)I + I*(l-a K (u+))X7(l-a K (u+))l} 



It is crucial that Z is a bounded and local observable. 

Theorem 16. II now follows immediately from (|6.2j) and Lemma 16.41 

The rest of this section is devoted to prove Lemma 16.41 We shall drop the labels ± on X and 
u etc for the rest of this section. 



6.1 General Properties of limiting states 

We now prove a number of results for the ergodic components of the limit points n. At this point 
T) depends on a macroscopic space point z, and a macroscopic time tj, and in principle also on the 
subsequence, but we will eventually see that n is in fact independent of the subsequence. 

Lemma 6.5 Let 7^,7 be normal states on a von Neuman algebra A, and 7„ — > 7 weakly. Suppose 
that A is a non-negative self-adjoint operator affiliated with A, such that "f n (A) is bounded by a 
constant M, uniformly in n. Then, hin^^oo 7 n (A) exists and satisfies 

7(A) < lim ln {A) 

n— ->oo 

Proof: Let A = j \dE\ be the spectral resolution of A. As A is affiliated with A, the projections 
Pk = J dE\ belong to A, and ~f n (A) = sup fc 7(^4/c), where A^ = AP^. The supremum is finite by 
the assumptions. Therefore, 

lim7„(A) = limsup7 n (A fc ) > suplim7 n (A fc ) = 7(A) 

n n k k n 

■ 

Recall that he and ne (we omit the superscripts +) are the average of the local conservative 
quantities Hi and N$ (|4.8|) . Let e and p be determined by 

e = lim uj(L*hiI), and p = lim uo{L*niI) (6-6) 

£^00 £— >oo 

Where necessary, we will indicate the dependence on uj by e(u>), and p(oj). In the following lemma 
we prove the existence and finiteness of these limits when the parameter a is fixed. 
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Lemma 6.6 For p— almost all states u, the limits e and p of (|6.6|) are finite. For the state 
n = j ' ui p(du>), we have the following bounds 

limsup 7/(7* hpl) < Ca~ 3 eo, limsup n(I*n£l) < Ca~ 3 po, (6-7) 



Proof: We have Av Av 7e "(/if) = eo, where eo is the initial total energy. This is a direct 

all boxes Q 

consequence of the fact that the energy is conserved by the dynamics. Therefore, for each box Q, 
we have 



Av 7e Qj (/i,) <Ca- 3 e 



By Lemma 16 . 51 with A = I*h£l, it follows that 



V (I*h e I) < lim Ay 7£ Qj (hi) < ]T Ay ^ (h t ) < Ca~ 3 e , 

which implies the bound for the energy (|6.7|) . The proof for the particle density is the same. As 
i] = f uj p(du>), it then also follows that e(ui) and p(uj) are finite for p— almost all uj. | 

Let A be a bounded observable in the local algebra .Aa C «4r3. E.g., A = j dx dyf(x,y)a+a y , 
where f(x,y) = unless x,y £ Aq. For concreteness, we assume that Ao contains the origin. We 
will also use the notation 

A ^ A)= w\L dxTx{A) 

Lemma 6.7 Suppose lim^^ = 0. For every translation invariant ergodic state u on A^3, any 
bounded local observable A and any continuous function f , we have the limit 

lunuj(r i f(AvA)I e ) = f(io(A)) (6.8) 

Proof: The proof rests on the following property of /: for A G -4a , we have 

I*t x (A)I = t x (A), ifr x (Ao)cA7 (6.9) 

Denote by A^ nt = {x £ A e \ t x (A ) C Aj}. Note that 

IMA^ nt | diamAo 

< 2% H =: <5 £ (6.10) 



|A,| " " I 

and that limg = 0. 

First, consider the function fix) = x. Then, using the property ()6.9|) . we have 



u{It Av(A)I) = f Tx{A )) + -L / r x (A) 

A f l A ^l JAi, nt l A d JA e \A^ 



I) 
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Without loss of generality we may assume uj(A) = 0. Using the definition of 5i Q6.1U[) . and the 
isometry property of I, we find 



Av(A)I) 

A f 



< (1 " k) 



1 



IA 



int. 



dxu(r x (A)) 



+ Si\\A\\ 



The two terms in the RHS tend zero, the first due to the ergodicity of u, the second because 6j> — > 0. 

Next, we prove by induction the result for fix) = x n , for all n > 1. Suppose we have the result 
for f(x) = x n_1 , i.e., 

limw(I*(AvA) n-1 I) = 

l A e 

Then, by the same arguments as above, we have the estimate 



lim 

I 



u{I*{kv{A)Yl) 

A f 



< (1 - Se T 


^((Av(A)f) 




Amt 



+ C(l-(l-fc)») 



and the result follows by the ergodicity of to. For arbitrary continuous functions /, (|6.8|) can now be 
obtained by approximating / by polynomials, uniformly on [— 1|^4||, ||j4||]. This proves the Lemma. 



Lemma 16.71 can trivially be extended as follows: 

Corollary 6.8 For any bounded local observables X,Y,A £ Aa > o,nd continuous functions f and 
g, we have that 



lim 

i 



u(PXf(AvA)Yg(AvA)I) - f(u(A))g(u(A))u(I*XYI) 

A e A e 







6.2 Extension to unbounded conserved quantities 

Lemmas 16.71 and Corollary 16.81 are general properties of ergodic states applied to bounded observ- 
ables. We now show how the one-phase region cut-off functions, which depend on unbounded but 
conserved quantities, can be included. This is a difficult step and we will have to use the spe- 
cial forms of the conserved quantities. The key technical estimate is contained in Lemma 16.101 
We remark that a naive application of Schwarz' inequality to prove Lemma 16.91 would produce 
expressions with six or more creation or annihilation operators about which we have no control. 

Lemma 6.9 Let r\ be any limiting point of {7^ J | e > 0}, let X be one of the components of w, 
and let Xi the averaged version of X. Then the following limits vanish: 



lim 77 (PB e X e [o-Uhi)a%(ne) - a^(e)a^p)]I) = 

'^00 

lim V (PB e X e [o-Z(n e ) - <r$(p)]I) = 



(6.11) 
(6.12) 
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for Bi = 1 or 

B e = a^(h e )a£(n £ ) 
Here e = e(rj) = J e(u>) dp,{uj), and similarly for p. In particular, we have 

lim V (r\a K (h e ,n e )X e a K (h e ,n e )} - \a K (e, p)X l a K (e, p) 1 i) =0 (6.13) 

and the same result holds if a K (hi,n£)Xia K (hi,ni) is replaced by a K (hi,ni)Xi or by Xia K (hi,ni). 

Proof: We start with the case Bi = 1. 

Recall, o K {e, p) = eri(e)<T2/j). There exist bounded functions <rf and such that erf (x)— af (y) = 
(x — y)af(x, y), for i = 1,2. Using these functions we can write 

a?(h e )o-%(n e ) - o-1{e)a^{p) = (af(h e ) - a^{e))a^{n t ) + o%{e){p%[n t ) - a%{p)) 

= o-%(n e )a1(h e ,e){he - e) + aZ(e)a$(n e , p){n e - p) 

Therefore, for a suitable bounded function /, for any ergodic state uj, we can write 

Lu(I*o- K (h e ,n e )X e o- K (h e ,n £ )) = Xa%(N e )f(h e , e)(h e - e)I) 

= u(lXaS{nt)f(ht)[hf-e B ]l) 

+u(l*Xo$(nt)f(ht)(h e -hf)I 



(6.14) 



where 



1 



I A, 



B 



with HX t = H Ai Ind(H Ae <B\A e \),so that \\H£ t \\ < B\A e \, and H" t ] H Ae as B oo. Introduce 

and use Schwarz' inequality to obtain 

Lo(l*Xo$(n t )f(h t ) [hf-e] l) 

< u[l*Xa^n e )f(h e ) [hf - e B ) i) | + (e B (u) - e(a/)) u(P<r%(n e )f(h £ )I) 

< 6u(i*x KM) 2 f(n e ) 2 x*i) + r^(r [/*f -e B ] 2 /) 

+(e B (o;) - e(uj)) ui(Po-%(n e )f(h e )I) 

Now, we integrate over uj with respect to the measure p, and take absolute values. The first 
term is uniformly bounded in £. As hf is bounded, the integrand of the second term vanishes 
for each uj, in the limit t — ► oo, by Lemma 16.71 As the integrand is bounded uniformly in uj, the 
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integral vanishes as well. The third term we use the argument of (|6.6|) to show that it vanishes in 
the limit B — >• oo: 

limsup \rj(hf - h e )\ < lim sup limAv 7£ Q ' j (/if - h t ) < Ca' 3 lime d \-/ £ t (H B - H)\ 

B B £ Q £ 

The RHS is independent of t and e, and vanishes as B — > oo. 

For the second term of Q6.14JI . we first apply Schwarz' inequaltity: 

u(l*Xa$(nz)f(h e ) [hi-hf] i) 
< &j(rx<r5(n t )f(h t ) [hf - h t ] f(h e )a%(n £ )XI^ + r 1 ^/* [hf - h t ] /) 

As before, the last term vanishes in the limit B — > oo. Since (he — hf) < hi, the first term is 
bounded by 

u(l*Xa2(n e )f(h e )h e f(h e )a%(n e )XI^ 

As / is bounded, we have that f(he)h(f(h() < Chg. Therefore, after integration over u>, and with 
the use of Lemma I6.1UI we obtain the bound 

5(Cri(rX{hi+ni)I) + C) 

which can be shown to be bounded in terms of the corresponding expectation in 7^, as before. 

In conclusion, as B and 5 are arbitrary, we have proved (|6,11|) for Bn = 1. It is straightforward 
to adapt the argument to prove also (|6.12l) and the case Bi = <rf (h^a^ (ni)- | 

6.3 Basic Estimate 

In the previous proof the following lemma was used. It provides a bound on the Hamiltonian 
sandwiched by bounded operators. 

Lemma 6.10 For fi- almost all translation invariant ergodic states oj, and Xp the averaged version 
of one of the components of w M (which are all self-adjoint), we have 

lo (rXta$(ne)hta5(ne)X t I) < Clu (I*[h e + n e ]I) + C (6.15) 

where the constant is independent of e,£ but may depend on a,M. 

Proof: Since X^ is particle number preserving, X^ commutes with ri£. Therefore, we can rewrite 
rewrite the quantity we need to estimate as 

lo (rX e a^(n £ )h £ a^(n £ )X e I)) = to (I*a^(n £ )X e h e X e a^(n £ )I) 
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hi is the sum of two terms, a kinetic energy and a potential energy term, which we wil treat 
separately. 

First, we consider the kinetic energy term: Va+Va^, defined with periodic boundary con- 
ditions. We start from the identity 

X e VatVa x X e = Va+X e X e Va x + Va+X e [Va x , X t ] + [X t , Va+]Va x X e (6.16) 

Note that Xg is a linear combination of linear and quadratic terms in M a v ,M ( see 1)2.613.513.80 ). 
Therefore, commutators of the form [a x M a yM ■> a x\i i a t M a y,M, Va+ ], etc., are bounded operators. 
More precisely, there is a constant Cm- such that 

II [X t , at] || < C M r 3 , and || [X£, Va+] || < C M r 3 . (6.17) 

These bounds will be used repeatedly in the following estimates. E.g., applied to the first term of 
(|6~TK|). they yield 

Va+XXVa x < C M Va+Va x . 
To bound the second and third term we first apply Schwarz' inequality: 

u^I*a^(nt)[X e ,Va+]Va x X e a^(ne)l) 

< 5Lu(l*a%(n £ )[X £ , Va+][X t , Va+]V 2 K (n^)/) + ^(/^(n^Va+Va^a^)/) 

The first term of the RHS is bounded and the last term can be re-absorbed into the quantity we 
started out the estimate. Thus, for the kinetic energy term and any of the Xt, we have an estimate 
of the form 

u(l*<T§(ne)Xiho^X t (^(nt)l\ < Cu}(l*cr$(nt)ho,t(T$(nt)l) + C . 
Similarly, for the potential energy we start from the identity 

Xl<l x (Xy dyClx 

= a x ayXeX e a y a x + a^a^Xt[a y a x , Xt] + [Xt, a x a y ]a y a x Xt 

and the bound 

a x ciy XgXta y a x < Ca x a y a y a x . 
For the commutator terms we have 

a x a y Xt[a y a x , Xt] = a x a y Xa y [a x , Xt] — a x a y X[a y , Xg]a x 
which can be estimated using Schwarz' inequality: 

2Reu(l*a%(nt)ata+Xtay[a x ,Xt]<T%(nt)l\ 

< u (i* a 2 (nt ) a+ a+ Xj a y a x a$ (n e ) ij + lo (i* a% (n e ) [a x ,X e ]a+ a y [a x , Xt] o% (nt ) 
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We < CMCtx a >y a y a x for the first term. The second term we use the identity 

[a x ,X]*a+a y [a x ,X} = a+[a x , X]*[a x , X]a y + a+[a x , X]*[a y , [a x ,X]] + [[a x , X]* , a+]a y [a x , X] 

The first term of the RHS is bounded by CMa y a y . The other two terms can be bounded by 
CMUydy + Cm by repeating the same procedure once more (first apply Schwarz' inequality, then 
use Ijfi.l7jl ). We conclude that 



dx I dyW(x-y)u(l*a$(n£)X£a+a+a y a x X£a$(ji e )I 

A e J A t ^ ' 

<C dx dy\W\{x-y)^u[ra^{ni)ata y v a y a x a^{nt)I S j 



A f J A 



+ u(l*a%(n e )a+a y a%(ne)l) + c} 
Now from the super-stability estimate, we have 



C I dx I dy\W\(x — y)a x a y a y a x < C I dx I dyW (x — y)a x a y a y a x + N \ 
J Ai J A e J A e J At 



Thus, 

r 3 J dx dyW(x-y) J dyu(l*aZ(n e )Xea+a+a v a x Xta$(n e )I^ 
Ci ~ 3 J a dx dyW(x-y) J dyu^I*a^(n e )a+a+a y a x a^n e )I^ 
+C dx dyW(x-y) / dyuj( I*a$ (ni)n^crf (n e )I 

JAp J At J ^ 



'At 

The last term is bounded. Combining these estimates, we have 



u[roZ(nt)X t htX t o%(nt)l) < Cu[l*a^(n e )h e ai(n e )Ij+C 
Since hp < h.£ + Cri£, hi + Cn^ > and [h^, ri(\ = 0, we have 

u(l**$(n t )h t *$(n t )l) < u(l*[h t + Cn e ]y 2 cTZ(n e )a%(ne)[h e + Cn,] 1/2 /) 



< Cu)\I*[hi + ne\I \ +C 
We can prove that uj(l*[he + ri(\l\ is bounded by using Lemma 16.51 



38 



6.4 Proof of main ergodic lemma 

We can now prove Lemma 16.41 

Proof: Recall the decomposition of r/ into its spatially ergodic components: 




Since X is bounded, by Lemma 16.91 there is £l Rt x(£) such that 

- a K )X e (l - a K )I) + u)(ry/a K (2 - a K )X y / a K (2 - a K )I) 
= (1 - a K {u)) 2 uj(I*X e I) + a K (to)(2 - a K {uj))uo{I* XI) + n K>x {£) 

where a K {ui) = cr K (limiu)(he),limguj(ni)). Therefore, 

7](rx e I) - r/(J*(l - a K )X E (l - a K )) - rj(I*^/a K (2 - a K )Xs/a K (2 - a K )I) 
= J n(dbj)[l - (1 - <J K {uj)f - a K {uj){2 - a K (uj))]uj(I*X e I) 

- J »{du)a K (u)(2 - a K (u))u{I*(X - X t )I) + Sl^X) 

As 1 — (1 — x) 2 — x{2 — x) =0, the first term vanishes identically. The middle term vanishes by 
the hypothesis that the only ergodic states of finite specific relative entropy in the one-phase region 
are the Gibbs states. The support of the function a K {uj){2 — (j k )(oj)) is such that only these Gibbs 
states contribute to the integral. The integrand vanishes by the definition of X (|4.1|) . since we have 
u>(X) = X(lim.gu(ug)). This concludes Lemma 16.41 | 



7 Relative entropy estimate 



We now summarize the estimates on the relative entropy we have so far. For any < T < Tq, we 
write 



s(7t I w !)| t - £ -i r = £ lT Av ^tA^ I w *< 

' L — t J n^4^. — 1 T< rlt 



d 

O^Ke-'T dt 

We compute the rate of change of entropy by (j3.4j) , (|3.10|) and 1)3. 16|) to have 



s(lt I cjf)| t=£ _ 1T = T q< Av ir {Tr 7t G(A £ ,a+ ,a M ) - e 2 d t log c £ (t)} + £ l M 

where G is defined in (EHUl) and £ l M < Ce~ cM ^ (l3~TnT) . 

Recall the meaning of the various length scales and cut-off parameters: e is the ratio of the 
macroscopic to microscopic length scale, M is the high-momentum cut-off, t is the length scale 
in the isometry I employed to define commuting local versions of the conserved quantities, a is 
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a length scale for averaging needed to make use of local ergodicity, k is the length scale used to 
smooth the characteristic function of the one-phase region, and 5 is a small parameter used in 
applications of the entropy inequality. 
Recall the convention 



3 3 3 

3 

3=0 fi=0 j=0 



We now apply Theorem 16. II to estimate Tr 7 t G(X £ , a~j^, om) by 

Av Tr 7t G(X E , at,a M ) <T X + T 2 

t<T/e 

where T\ and T 2 are defined as follows: 

Ti = - Av 7t Av ( VX(et,ex)» I*(a K w Mx a K )l 

0<t<e~ 1 T x \ ' ' 

T 2 = - Av lt Av(v\(et,Ex)*I*{{l-a K )(u x>e )w Mx {l-a K )(u x>e )}l 

0<t<e- 1 T x \ 

We need to compute w which we state as the following lemma. 

Lemma 7.1 We have the following identities 

w^ = A^ j = (),•• -3, M = 0,--- ,4. 

where the functions Aj are given in (|2.15|) . 

These relations follow directly from the definition of the Gibbs states, the expressions for Aj in 
(|2. 15|) , the calculation of the currents in Section El and the virial theorem proved in Section I1UI 
By construction of the wf, the time derivative of logc £ (i) can be expressed as 



d f 
—e 2 \ogc £ {t) =e 3 / dx{d t \ ■ q)(et,ex) 



where q is the solution of the Euler equations that we are considering. Recall the following identity 
about the Euler equations: 



Recall also Aq(ci) = q. We can rewrite 



3 

J J2 A M*)) •V i A(q(X))dX = 0. 



d_ 



e 2 log c e (t) = e 3 J dx(d t X ■ q)(et,ex) = e 3 J dx(VX ■ A)(et, ex) 
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Together with Lemma 17. 11 we have 
Tt - i 2 logc £ (t) = - Av lt Av \v\(et,ex)» { I* (d K ± M x d K ){u x/ )I - w(q(et,es)) } 

(it 0<t<e _1 T a: L ' 

Denote 

llVAHoo = ||VA||oo + II^AHoo 

and introduce the functions 

Im(A, u) = VA. [ (a K w M a K )(u) - w(q) ] , 

T 2 (A,u) = ||VA|| 00 [(l- CT K )(u)(/ i + n)(l-a ft )(u)] , 

where h, n are the energy and density components to u and q is the dual variable of A defined in 
iprTTj) . 

We can bound w M in T2 by the cutoff Lemma 15. II Thus we have 

Ti + Ta-^e 2 logout) < Av lt Av fl* {-T l M + MT 2 }(X(et,ex), u x/ )l) 

dt 0<t<e~ 1 T os \ J 

Therefore, we have 

s{lt\u e t )\ t , T <T Av 7t Av(r{-T 1 M + MT 2 }(X(et,ex),u x/ )l)+£ 1 M (7.1) 
where £ l M < Ce~ cM2 (Hmi) . 

7.1 Reduction to large deviation 

Recall the standard thermodynamics pressure is defined by 

VKA) = lim r 3 io g z, iA 

£^00 

Define the entropy 

s(q') = sup[Aq - ^(A)] 
A 

and the rate function (notice we also use I for the embedding into the standard torus A e -i) 

/(q / ,A) = S (q / ) + ^A)-A-q / 
The rate function has the following property 

/(q',A)>0, /(q,A)=0 
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where q = dip(X)/d\. Furthermore, if the Gibbs state with chemical potential A is in the one 
phase region, we have 

Hess J(q, A ) > ct 

for some c > 0. 

The main large deviation estimate we shall use is given in the following lemma. This lemma 
will be proved in Section \tt.dl 

Lemma 7.2 Suppose X is a bounded smooth function so the Gibbs state with chemical potential 
X(x) is in the one phase region for all x. For any bound smooth function G that satisfies the 
condition 

|G(A, q)| <C{e + p) (7.2) 

where e is the energy and p is the density. Then there is a 5q > depending only on C and a 
convex functional I such that for all < 5 < 5o 

I(q,A)=I(q,A) 

in a small neighborhood o/q = dip(X)/dX and 

lim lim 7 Av (i* G(X(ex),u+ )I 

< [ dX sup [G(X(X)), q \X)) -5- 1 I(q / (X),X(X))] + r 1 lim S (7l^i) 
Here the sup is over all functions q' (X) . 

7.2 Conclusion of the relative entropy estimate and proof of the main theorem 

We now apply Lemma 17.21 to estimate (j3.17|) . Since we need the bound (|7.2|) . we set G = 
M~ l {-T X M + Mr 2 }. Thus we have for any 5 < 5 



where 



- Av 7t Av (P{-r L M + MT 2 }(X(£t,ex),u+ e )l) < i? 6 + 5- l Ms(j t \ ojf) 

0<t<e~ 1 T x V ' / 

R 6 = [dx sup {{-Tl[ + MT 2 }(\(x)),q[(x)) - r J M J(q'(a;), X(x) ) ) (7.3) 

J q' 

where / is related to the rate function defined in Lemma 17.21 We now estimate the dependence of 
rjj, on M. 
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Lemma 7.3 There is a constant c > such that 

ri f (A,q')=r 1 (A,q')+e- cM2 

where 

rH\ q' ) = VA. [ ( ff > K )(q') - w(q) ] 

This lemma can be proved following the idea of the proofs of Lemmas 13.21 and 13.31 It is part of 
our assumptions that the Gibbs states satisfy the cutoff assumptions. 

Now, we can conclude the relative entropy estimate and the proof of Theorem 12.11 
Proof of Theorem I2TT1 

Recall q = dip(X)/d\. Clearly, r 1 (q(et, ex), X(et, ex)) = 0. The first derivative 

dT 1 (X(et, ex), q(ei, ex) ) _^ 
<9q(ei, ex) 

is equivalent to the Euler equation as checked in ^Sj. Recall T 2 (A(X) , q'(X) ) is nonzero only 
when q'(AT) is away from q(A"). Thus we have for |q'| < C 

{.r 1 + Mr 2 }(A(X) , q'(X) ) < CM(q'(I) - q(X)) 2 

Furthermore, from the definition of Tj we have for all q' 

{-r 1 +Mr 2 }(A(X),q / (X)) < CM( |q'|(X) + 1) 

Since I(q.'(X) , \(X) ) > and I(q.'(X), \(X) ) = only when q'(X) = q(X), for 5 small enough 
we have 

sup [{-r 1 + Mr 2 }(q'(l) , X(X) ) - 5- l MI(c({X) , \(X) ) ] < e~ cM2 
q' 

We thus have 

s{lM) < S~ x Me f 
Jo 

By intergrating this inequality is (i.e, using Gronwall's inequality), and using the fact that et < Tq, 
we arrive at the bound 

s^M) < 5- l MT,e & - lMT ° 

Taking the limits lim K ^o lim tt ^o hni^oo lim e _>0) we get the inequality 

lims( 74 |u; t £ ) < C8- l MT e ^MT -cM\ 

We can now let M — > oo and conclude the proof of Theorem 12.11 We emphasize that we need the 
error term stemming from the high-momentum cutoff to be smaller than e~ CM for any C > in 
order to have our results hold for t < CTq for arbitrary To. This is guaranteed by the Maxwellian 
bound in the cutoff assumption II. 1, expressed by Q2.19[) . 



t>\<A', 



+ Ce~ cM2 + n M (£, k) + n{e, a) + C M e 



d£ 



Ce~ cM + n M (£, k) + fi(e, a) + C M e~ a 
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8 Thermodynamics and large deviation 



We now prepare the way for the proof of Lemma l7.21 Our approach to large deviations for quantum 
Gibbs states and local Gibbs states is quite different from the explicit analysis in for the ideal 
gases. We first introduce the following local Gibbs state with independent subcubes. 

8.1 Local Gibbs state with independent subcubes 

Divide the torus A = A e -i into unions of non-overlapping cubes of size i. To fix the grid, we assume 
that the origin is the center of one small cube. Denote a typical cube by a. Recall the configuration 
space S(a + ) and define the configuration space 

S(A«) = ® ae ffl3nA5(a + ) 

An element in this configuration space can be denoted by 

x = (• • ■ , (oij , Xj), • • • ), Xj 6 o.j 

The (Fock) function space r(A<+)) is the I? space of antisymmetric functions on 5(A^ + ^). Notice 
that 5(A(+)) + S(A+) and r(A«) ^ T(A+). 

Recall A = A £ -i is a torus. Define l£ from T(A) to r(A^) (cf @]) by 



The crucial fact is that 1 is an isometry. 
Lemma 8.1 X is an isometric embedding, i.e., 



ip(x) . 



Proof: Recall from the construction of x the- relation (|4.2|) implies that 

aeH 3 

We can prove the isometry by the following identity: For any two wave functions / and g, we have 



a 3 



<foj|Xcy(a?j)| 2 



/(xjg^) = (/, 5 ) 



For isometric embeddings, we have the following useful bound. 
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Lemma 8.2 Suppose X : Hi — > TI2 is an isometric embedding. Then 

Tr Hl e M < Tr n 2 e A 

Proof: We can assume that H\ is just a subspace of H2 and X is the natural embedding. Let 
be the othornormal eigenvectors of A in 7i\. The the claim follows from the following Peierls' 
inequality: Suppose (f>j are othornormal. Then 

^ e i*iMi) < Tr e A 



The following Lemma shows that I* XI = X*XX for a suitable class of observables. 
Proposition 8.3 Suppose X is the observable 

X = J dx 1 - ■ ■ dx k dyi - ■ ■ dy k f(xi, . . . ,x k ;yi,. . . ,y k )a+ ia - ■ ■ a+ ka a ykiC( - ■ ■ a yua , 
Then we have I* XI = X*XX. 

Proof: The following identity is a direct consequence of the definition of X: For n > k, 

j Zn ) 

= S ai>a --- 5 akt aX a (yi) ■ ■ ■ Xa(yk)Xa k+1 (Zfc+l) ' ' ' Xa n 0„) 

It follows that, for any (f>, ip E L 2 (A xra ), n > k, 
(4>,X*XIip) = / dxi---dx k dyi---dy k f(x 1 ,...,x k ;yi,...,y k ) 

Ja+ 

= ^ dzk+i ■ ■ ■ dz n dxi ■ ■ ■ dx k dyi ■ ■ ■ dy k S ai>a ■ ■ ■ S akia 

ai,...,a„ 

XXQ i+ i(%l) 2 • • • Xa n (z n ) 2 Xa(xi) ' ' ' Xa(xk)Xa{yi) ' ■ ■ XaiVk) 



xf(xi, ...,Xk\yi,-. .,yk)(/>(xi, . . .,x k ,z k +i, ■ ■ .,z n )if>(yi, . . .,y k ,z k+ i, ...,z n ) 

The sum over a\, . . . ,a n can be carried out using the Kronecker delta's and l|8.1|) . As n > k, (f>, 
and if) are arbitrary, we have I* XI = X*XI by the formula of / in (|4.4|) . | 

We now construct a "special local Gibbs state". Recall that u a + is defined in (|4.8|) and its 
component commute. For a smooth A, let tDt' be the state 

£o £ y = Tr — — exp e~ 3 Av X(ea) ■ I*u a + I 
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where the average of a is over a £ 11? Pi A and c £ g(X) is the partition function defined by 

c £ i{X) = Tr exp e~ 3 Av X(ea) ■ I*u a + I 

L a 

Here the trace is over T(A^). 

Assume for the moment we can drop the / and the small cubes are independent. Then c e /(X) 
can be computed easily. The following Lemma asserts that this is essentially correct. Recall the 
partition function defined in (|2.1U|) . 

Lemma 8.4 

lim lim (e 3 log c £ (X) - Aw T 3 log Z e (X(ea))) = lim lim (e 3 log c £ t (X) - Aw T 3 log ZAXiea))] =0 

(8.2) 

Proof: Upper bound: We first state an upper to e 3 logc e ^. Notice that it is an inequality with no 
limits or other constants. 

e 3 logc M (A) < Avr 3 log^(A(ea)) (8.3) 

a 

From Lemma 18.31 we have 



logc £i £(A) = logTr exp 



s~ 3 T Av X(ea) • u„+ 1 



If we can neglect I* and I, then ()8.3|) follows from the fact that different cubes are considered 
independent. Lemma 18.21 shows that we can remove I* and I to have an upper bound. This 
concludes the proof of (|8.3|) . 

Lower bound: Since A is fixed, we shall drop it in the subscript. Consider the entropy 

< s ( uj £ | uj £ ' e ) = R 2 + e 3 log c £ ,e - e 3 log c £ 

where 



R2=io,(e 3 



12 = iv £ yt ' J dxX(ex) ■ u x — Av A(ea) • I*u a + I) 
From Lemma 14. 11 we have 



We have thus proved that 



lim lim R 2 = 



lim lim [ e 3 log c £ti - e 3 log c £ ] > (8.4) 



To conclude Lemma 18.41 we now obtain a lower bound on c £ . This is the standard procedure on 
the thermodynamics and we shall give only sketch. We first divide the cube of size e~ l into cubes 
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of size £(1 — VI) with corridors of size 2y/I. Now we impose the Dirichlet boundary conditions on 
the boundary to obtain au upper bound on the kinetic energy. The partition function is bounded 
below by restricting the configurations so that there is no particle on the corridors. Now there is 
no interactions between different cubes and we obtain a lower bound of c e in terms of average over 
Dirichlet boundary conditioned partition functions in cubes of size i(i - VI). Since we can take 
77 — > after £ — > 00 and partition functions is independent of boundary conditions, we have thus 
proved that 



lim lim^logCe > Av £ d log Z e (X(ea)) 

£— »csoe— >Q a 



This concludes the Lemma. 



8.2 Large deviation for commuting variables 

Recall that u a + is defined in Q4.8|) and its component commute. We shall take x = and de- 
note u a + by u^. The following Lemma is a standard application of large deviation theory (or 
thermodynamics) to commuting variables. 

Lemma 8.5 Suppose X is a fixed constant so the Gibbs state with chemical potential X is in the 
one phase region. For any bound smooth function G satisfies that 

\G(X, q)|<C(e + p) 

where e is the energy and n is the density. Let be the finite volume Gibbs state defined by 

u\/(X) = ) Tr exp [ (A, u^) A J X 

with periodic boundary condition. Since the components of U£ are commuting, we have 
1 



-Tr exp [(A, u e ) Ae + 5£ 3 G(X, u e )] = <jj\£ ( exp [8£ 3 G(X, u e 



Then there is a 60 > depending only on C and a convex functional I such that for all < 5 < 5o 

I(q / ,A)=/(q / ,A) 
in a small neighborhood o/q = dtfj(X)/dX and 

lA^logwj^exp^GCA, u,)] ) < sup [SG(X, q' ) -I(q',A)] 

Here the sup is over all constants q'. 
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We first sketch the idea of the proof for Lemma 18.51 The rate function / can be understood in 
the following way: The probability to find the with a given value q' is given by exp[— |A^|s(q')] 
with the entropy given by 

S (q / )= S up[A-q / -^A)] 
A 

We now write Tr z \\) ex P [ U ) A <J as 



j dq'exp[|A £ |{A. q'-^A)-^')}] 



This gives the last variational formula. 

Proof of Lemma 18.51 We shall drop the constant parameter A in G in this proof. Since the 
components of u commute, we can define the joint distribution i^(du) of u w.r.t. the state co\ t £. 
Thus 

w v (exp [5£ 3 G(X, u e )) ) = J d^(u) (exp [^ 3 G(A, u e )] ) 
We now approximate the integral by the summation so that 

I dm(u) exp [5£ 3 G(\, u)] < P\A\ u ~ £m \ ^ £ ] ex P [8PG(em)] 



where 



G £ (y)= sup G(x) 
|x-y|<£ 



and P\ t £ denotes the probability of the event described in its argument, with respect to the state 
0J\£. We can bound the summation by 



e~ 5 j dx PxA\ u ~ x l < £ ] ex P [^ 3 G e (x) 



We have 

P v (|u-x| <e) < Px,e(Z - 

for all Notice that from the Chebeshev inequality we 
have 

iVU-u>£-x-|4k) <e-***+*®°J %(u) 
Let ipt{X) be the pressure defined by 

^(A) =r 3 logTr (e^ A ' u ) 

so that 



j dm{u) e^ u = exp {£ 3 [M£ + A) - AW] } 
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Thus 



for all 



iVU ■ u > £ ■ x - |C|er) < exp { - £ 3 [£ • x - + A) - e|£| + ^(A) 



A 4 + £ 4 > . 



In particular, we have 

P V (|u-x| <e) < exp{ -t 



i7- 1 <«j+A J -<7j- 1 ,j=0-3 



(£-x- ^(£ + A)- £ |£|+^(A))} 



sup 

A 3 <, 
?7<^4 + A4<T7 1 

The existence of thermodynamics states that 

lim U(A)-V(A)| =0 

£— >oo 

uniformly in compact interval away from A4 = 0. Fix a small constant r/ > 0. Define 

s v (x) = sup U-x-V(O) 

-rj- 1 ^-^!?- 1 , j=0-3 
»7<54<"7 _1 



We have 



where 



Define 



sup (£. x _ ^(£ + A)-e|£|) 

-'7- 1 <« 3 +A J < J? -l,j=0...3 
»7<«4+A4<'7 _1 

< S,(x)- A • x - C(?7, A)e + C £ 



lim Q = 

£^00 



I v {x) = s^(x) - A x + */>(A) 



Thus we have 



P V (|u-x| <e) <exp{-£ 3 (l r) (x)-C e -C( V ,\)e)} 
We now have the estimate 

r 3 log J dfi t (u) exp [5£ 3 G(X, u)] < r 3 log {e~ 5 J dx P x/ ( |u - x| < e) exp [Sl 3 G £ (x 



< C 3 lo 



■g j [dx exp [-f 3 (/ r ,(x) - <JG e (x))] +Q + C(r?, A)e + r 3 |loge| 



The error vanishes in the limit lim e ^o lim£-+oo- The integration can be calculated using the Laplace 
method to give 

lim r 3 log / dx exp \-£ 3 (l v (x) - <5G e (x))l < sup [5G e (\, q' ) -I v (q',X)] 

£->oo J L J q / 
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Clearly, we have 

lim sup [5G e (X, q' ) - I,(q', A ) ] = sup [SG(X, q' ) - j^q', A ) ] 

e >U q/ q/ 

We now collect some properties for and I v . Notice that, by definitions, and I v are still 
convex. Furthermore, we can check that if rj is small then 

/„(x) = /(x) 



in a small neighborhood of q. This proves the Lemma. 



8.3 Proof of Lemma EH 

Recall a indices disjoint subcubes of width £. By the entropy inequality (|3.17|) . we have 
- Av 7e Av (i* G(X(sa),u a +)l) <R + ^ 1 Ms( 7£ I Cb £ /) 

0<t<e~ 1 T a \ J 

where, for any 5 > 0, 



R = e d 5- l \ozTi 



exp 



e 3 Avl*< \(ea) ■ u a + - SG(X(ea), u a+ ) \I 



(8.5) 



c £ /(A) 

We first control the last term s{^ e \ u)l' f ~) by writing it as 

s( 7e | uf) = s{ le | + e 3 logc e (A) - £ 3 logc £/ (X) + R 4 

where 

i?4 = 7 e ( e 3 J dxX(ex) ■ u x — Av A(ea) • J*u a + /) 
From the Lemma l4.1( we have 

lim lim R4 = 

We now estimate R. Using the argument in the proof of Lemma 18.31 we can drop the operator 
I in (|8.5|) to have an upper bound. Since the cubes indexed by a are independent, we have 



R < <T 1 Av Q a - ( e 6 log c si (X) - Av £~ 3 log Z t (X(ea)) 



where 



Q a = £~ 3 log Tr ^(A(ea))- 1 exp 



£ 3 \ X(ea) ■ u a+ - 5G(X(ea) 



The last term ^e 3 logc e .^(A) — Av^ 3 log Ze(X(ea))J vanishes by Lemma 18.41 Notice that the 
components of u a + commute. Thus the trace is over functional of commuting operators and we 
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are essentially the same as in the classical theory. Thus we can apply Lemma 18.51 to estimate Q a . 
Summarizing, we have 

lim lim7Av ( I* G(X(ea),u a+ )l) 

< I dX sup \G(X(X)),q'(X)) -5- 1 I(q!(X),X(X))] + 5~ l Kmsfa I u\) +R 4 

J q'(X) 

Notice that the right side of the inequality is independent of the location of the grid. If we average 
the grid over the cube of size £, we can replace the left side of the inequality from averaging over 
a to averaging over all points x on the torus. This proves Lemma 17.21 

We now state a corollary to Lemma \7~7\ 

Corollary 8.6 Suppose X is a bounded smooth function so the Gibbs state with chemical potential 
X(X) is in the one phase region for all X . Suppose j e is a sequence of states such that the specific 
entropy s(7 £ | u x ) satisfies 

lims(7 £ | u x ) = 

e^O 

For any bound smooth function J on the unit torus, we have 

lim 7 £ e 3 J dx J (ex) ■ u(ex) = J dX J(X) ■ q(X) 
Proof: Since J is bounded smooth, from Lemma 14. II we have 

limlim 7 ( £ 3/' (fcJM .u M -AvJM.<,)=0 

i^oo e^o v J x x ' 1 ' 

We now apply Lemma 17.21 to have 

lim lim 7 e Av J(ex) ■ ( ut c — q(e^) ) 

^oce^O x v ' v x ' 1 v ' ' 

< [ dX sup \j(X) • ( q'(X) - q(X) ) - S^Ifa^X), X(X)) 1 + (T 1 lim s( le \u; x ) 

J q'(X) L J £^0 

Since / > and I(q'(X), X(X)) = only when q'(X) = q(X), the sup is bounded by C5. To see 
this, consider the model problem 

supx — 5~ 1 x 2 < 5 . 

X 

Recall the assumption lim^o sCTeI^a) = 0. Since we can choose 5 arbitrarily small, we prove the 
corollary. | 
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9 Calculation of the currents 



The current density operators w kx are implicitly denned by ()3.9|) . i.e., for any test functions J 
( J M ), fj, = 0, . . . , 4, they should satisfy 

3 terms containing second and higher 

> w fc,x)) — derivatives of the integrated with 
fc=1 densities of bounded expectation. 

If we apply the same sign convention for dot products with w as the convention adopted in 
for u, this means we are looking for the definition of w kx , such that, for any test function J, the 
following formal identity holds: 

i \ dxJ(x)[H, u%] = ^2 / dxS7 kJ{ x ) w kx + integrals with higher derivatives of J , (9.1) 
J k=i J 

where H is the formal Hamiltonian H = J dxh x , with h x as defined in (|2.6|) . 

In order to compute the commutators we use the canonical anticommutation relations (|2.4|) and 
integration by parts. The commutation relations involving derivatives such as V k a x etc., are most 
easily derived by taking derivatives of the appropriate commutation relations without derivatives. 
E.g., the identity 



8{x — v)a^a y — 5(y — u)a x a v 



follows directly form ()2.4|) and, by taking derivatives with respect to u, also leads to 

Vfca^a„, a x a y = 8{x — v)V} t o^ a y + $k{y ~ u)a x a v 

where 8k is the derivative of the delta distribution with respect to the fcth component. It is 
straightforward to derive all other necessary relations in the same way. E.g., 

[VfcO+VfcOy, a+a x ] = -8 k (x - y)V k a+a x + 8 k {x - y)V k a y a+ . (9.2) 

There are essentially three cases to consider: i) \i = 0, ii) \i = 1, 2, 3, and iii) \i = 4. 

i) fi = 0: calculation of : 
As n x commutes with the potential part of the Hamiltonian we only have to consider the kinetic 
energy term, which can be computed using ()9.2|) . After integrating by parts, we get 

i J ' dxdy J(x)[-VayVa y ,a x a x ] = J dxV J(x)-i[V k a x a x — a x Va x ] 

By comparing this result and we find agreement with the definition of w k as given in (|3.5j) . 
Note that, in this case, no higher order derivatives of J appear. 
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ii) n = j = 1,2, 3: calculation of wi . 

Now, both the kinetic energy term and potential energy term both yield non-trivial contributions. 
First, we compute the kinetic energy term. 

i j dxdyJ(x)[-VayVa y ,pi.} = - J dxdyJ(x)[Vja^a x , VdyVa y ] + h.c. 

where here and in the following h.c. stands for the adjoint of the preceding term(s). By using the 
commutation relations and integration by parts we find the following expression for this quantity: 

1 3 

dxdyJ(x)\pi, - ^Vfcd+ V 



2 ■ • / '"'' 

k=l 



— I dxJ(x) [Vja^AcLx + Aa x Vja x ] — V k J(x)Aa^a x + h.c. 

J k=l 

y~] / dxJ(x) [VfcVja+Vfca^ + V k a+V k Vja x ] 
k=i 

+2^ J dxV k J(x)Vja^V k a x -J dxVj,J(x)Aa+a x + h.c. 



k=l 

After further integration by parts and reorganization the result can be written as 

3 

L 2 



dxdyJ{x)[^ V k a+V k a y ,p> x ] 
1 k=i 

3 r i i 

= / \ / dxV k J(x)- [Vja^V k a x + V k a x ~V j a x ] + -V fc V,J(x) [V^a+a^ + a+VfcO^.] 
fc=l ^ 

The first term of the RHS in this expression determines the first term (|3.6|) . Note that this time 
higher derivative terms appear that are not included in the definition of wi , but they contribute 
to the error terms. 

To calculate the contribution from the potential energy term in the Hamiltonian, we start from 
the identity 

V ja^a u , a x a y a y a x = 5(x — u)V ja^a y a y a x + 5(y — u)a x \/ja^a y a x 

+Vj5(y — u)a x a y a u a x + Vj5(x — u)a x a y ~a y a u 
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1 



which leads to 

dxdyW(x-y) J duJ(u)-[V ja^a u — a^Vja u ],a x a y a y a x 

dxdyW(x — y) J(x)[Vja^a y a y a x + h.c] + J(y)[a x V ja y a y a x + h.c] 
+ [Vj J{x) + V jJ(y)]a^a y v a y a x 



(9.3) 



/ 



dxdy 



J(x)Vj, x W(x -y) + J(y)V^ y W{x - y) 



Due to the spherical symmetry of the potential we have 

VjW{x -y) = W'(x - y) 



(x - y)j 
\x - y\ 



(9.4) 



Using this identity we can write 1)9.3)1 in the form 



(J(x)-J(y))W'(x-y) 



- y)j 



(X x &y QjyQj X • 



F - v\ 

As the range of W is finite by assumption, we can Taylor expand J(x) — J{y) to rewrite this quantity 
in the following form: 

— / dxdyVkJ(x)W'(x — y) j j -a x a y a y a x + higher order derivatives of J . (9.5) 



fe=l 



f - y 



Recall that, by definition, only the coefficients of the first order derivatives of J are included in the 
w tensor. Therefore, combining ()9.2|) and ()9.5)1 and also including the appropriate factors 1 /2 and 
- signs, we find the expression for claimed in ()3.6j) . 

iii) = 4: calculation of x . 
The calculation of the energy current proceeds in the same way as the previous cases, but there 
are more terms and terms with higher derivatives. The contribution from the kinetic energy in the 
Hamiltonian to the kinetic energy current is, up to a trivial constant, given by 

i / dudx.J(u)[Va^Va u , Va x Va x ] = i / dudxJ(u)5k t i(u — x) Vfca^Vza^ — Via x Vkd u 

3 



k=l 



^2 / dxV k J(x) V fc a^Aa x — AaJVfca a 



where Si-i is shorthand for VfcV;<5. This yields the first term of the energy current. 

The potential energy term in the Hamiltonian does not contribute to the potential energy 
portion of the energy current due to the fact that the following commutators vanish: 



ci u a v d v d u ^ a x a y ci y (i x 
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To calculate its contribution to the kinetic energy current we start form 
J J{u) \Va+Va u ,a+a+a y a x 

3 f r 

= Y duJ ( u ) 5 k (u - x)V k a+a+a y a x + 5 k (u - y)a+V k a+a y a x 



k=l ' 

-S k (u - y)a x ayV k a u a x - 5 k (u - x))a x ciy a y V ' k a u 
= -V x (J(x)Va x )aya y a x - a x V y (J(y)Vay)a y a x 
+a x a y V y (J(y)Va y )a x + a x a y a y V x (J(x)Va x ) 

By multiplying this expression by W(x — y), and integrating over x and y, and integrating by 
parts, we find 



dxdyduW (x — y)J(u) Va„ Va M , a x a y a y a x 



[ J{x)V Kx W(x-y) 
k=i J 

3 

+ £_/ J(v)Vk, v W(x-y) 



V k a+a+a y a x - a+a+ayVkaa; 



atVka y a y a x - a+a+Vfca^ 



The contribution of the kinetic energy to the potential energy current is obtained in a similar way. 
The result is 



/ 



dxdyduW(x — y)J(x) a^aya y a x ,Va^Va v 



k=i 

3 



= / V k , x (J(x)W(x - y)) V k a+a+a y a x - a+a+a y V k a x 



atVka y a y a x - a^a^V k a y a x 



-Y, I J{x)Vk, y W{x-y) 
k=i J 

The last four terms become the last two terms of the energy current: 
i j dxJ(x)[h x ,H] 
i 3 f 



k=i 



i f 

~Y (V k J)(x)W(x-y) V k a^a^a y a x - a^a^a y V k a x 



k=i 

3 



( J ( y ) ~ J {x))V Ky W(x - y) a+V k a+a y a x - a+a+V k a y a x 

k=i J 
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By the same argument as for (|9.5[) . the last term can be rewritten in the form 

f V k J(x) \W(x - y) ( x -y)®( x ~y) ] \ a + Va + ayax _ a +a+Va y a x ] + 0(J" 

J 1 \ x y\ J 



-i 

k=l ' 



10 The virial Theorem 

The purpose of this section is to relate the expectation values of the RHS of the dynamical equations 
Q3.9|) in a Gibbs state with specified values of the densities of the conserved quantities (local 
equilibrium), to these quantities themselves in order to obtain a closed set of equations. To achieve 
this we will make use of canonical transformations relating Gibbs states with respect to reference 
frames with different velocities. This will allow us to use reflection symmetry of Gibbs states at 
zero total momentum. A second element we will need is the Virial Theorem to relate the so-called 
virial to the thermodynamic pressure. We start with the latter. For the convenience of the reader 
we first recall the main definitions. 

Consider a system of particles in a finite volume A C M. d , interacting via a pair potential W. 
The pressure at inverse temperature and chemical potential fj,, P{(3, /x), is defined by 

P(f3,u)= lim — i-logTr e -P{H ,A+V A -»N K ) ^ 
A^R d p\A\ 



where 



A 



2 

Va = ^ / j drdi) lT(.r - ■ i))apiJo ,,(!,- 
Na = / dxa^dx 



v 

AJA 



J A 

The trace is taken over the Fermion Fock space with one-particle space L 2 (A). For our purposes, 
we can simply consider A to be a cube of side L, and define the operators with periodic boundary 
conditions. We will write Va(W^) when we wish to indicate the pair potential function explicitly. By 
our general assumptions, the limit (|1U.1|) exists and we will restrict ourselves to the one-phase region 
of the phase diagram. In particular we assume that the pressure is continuously differ entiable. 

Gibbs states at non-vanishing total momentum are defined by introducing an additional La- 
grange multiplier for the momentum as follows. 

u x (X) = lim ^-Tr Xe^^^-"^-^ (10.2) 
A^R d Z(\) 
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where A = ((3,a,fi), a = (ax, a-i, 03), are constants, and Pa is the total momentum operator in 
the volume A defined by 



Pa = - dxVaJa x - ajVa x . 

2 -/A 



and 

Z(X) = Tr e ~l 3 ( H o,A +v ^~ aP ^~t lN ^ 

is the partition function. 

The kinetic energy density is defined by 

ekin(A oc,(Jt)= lim r^Up * JHq a) 
A^R d \A\ 

The limits A — > R 3 exist and are independent of the boundary conditions under general stability 
assumptions _18 ]. We will use the abreviations ujp^ = wp,o,n and e^i n (f3, fj,) = e^nCA 0, /i). 
The virial of the potential W in the volume A is denoted by Va(W) and is defined by 

Va(W) = \ I I dxd V VW(x -y)-(x- y)ata+a y a x (10.3) 
1 JaJa 

and the density of the local density virial is given by 

v x = - / dy VW(x - y) ■ (x - y)a x ~a+a y a x (10.4) 

As we have assumed that W has compact support, u x is well-defined. 

Due to Galileo invariance, the Gibbs states for different values of a are related by a canonical 
transformation, which is why in the statistical mechanics of global equilibrium situations the total 
momentum is usually assumed to vanish. 

The canonical transformations relating u>f) u and the states up a u, are defined as follows. Let 
s 6 M. d , and consider the unitary U s on L 2 (M. d ,dx) defined by (U a ip{x)) = e ls " x ^/>(x). The second 
quantization of U s implements an automorphism 7 S on the Fermion algebra given by 

7a (a(/)) = a(U s f), ls (a + (f)) = a + (U s f) . 

One can easily verify that the action of 7 S on the operator- valued distributions a x , Va x , and their 
adjoints, is given by: 



7s (a x ) = e JS ' x a x , 7s (a+) = e ls ^al 

7s (Va x ) = e~ is - x Va x - ise'^a^ 7 s (Va+) = e is ' x Va+ + ise is ^a 



X 



Clearly, 7s 1 = 7 _ s . 
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With these relations it is easy to check that 

7s(Va x Va x ) = Va x Va x + |s| 2 a x a x + is ■ (a x Va x — V x a x ) 
Hence, the kinetic energy transforms as follows: 



In the same way we see that 



It follows that 



7s (#a,o) = #a,o + ~|s| 2 ^Va - s P A (10.5) 



7s (iV A ) = N A (10.6) 

Js(Pa) = Pa-sN a (10.7) 

7a (^v) = V A (10.8) 

7*(Va(W0) = V A (W)) (10.9) 



7a ( e -^A-W) = e 



By putting s = a, replacing \x by [i — ||ct| 2 , we obtain 

e -/3(tfA-«-P A - M AT A ) = ^^-/3(P A -( M +I| a | 2 )JV A )^ 

As the trace is invariant under 7ct , this implies 

Z{\) = Tr e -P(H A -^P A -^N A ) = ^ e -/30f/ A -( M +iH 2 )7V A ) = Z 0^ 

where, for A = {(3, a, fx), we define A = (/3,0,/i + ||a| 2 ). Using this relation between partition 
functions and the invariance of the trace under canonical transformations, we immediately get 

"A(7«(*)) = ^TJ Tr J a (X)e-^- a - p ^ N ^ 

= _L T r x e -/3(ffA-(A*+§l°mA) = W3 .m (10.11) 
Z(A) A 

In combination with (jlU.6f) this implies 

u x (N A ) =u A (N A ), 
and hence p(A) = p(A), and we will simply write p. Also 

kj a (Pa) = aLu x (N A ) • 



58 



If we apply this to X = -Ha,o and combine this with Q1U.5[) . to relate the kinetic energy densities of 
lu\ and wt, we find 

e kin (/3,a,/i) = e kin (l3,p + ^|«| 2 ) + ^\^\ 2 p , (10.12) 

where we have also used tu^(P\) = 0. The relation IjlO.lOj) between partition functions immediately 
implies the following property of the pressure: 

P((3,a,p) = P((3,0,p+^\a\ 2 ) (10.13) 

One interpretation of this relation is that the chemical potentials at different values of a, when 
regarded as a function of the particle density p, satisfy 

Ha(fi) = Mo(p) - ^|"| 2 • 
We can now prove the virial theorem in the form we need. 

Theorem 10.1 (Virial Theorem) For a three-dimensional translation innvariant system with a 
continuously differentiable pressure function, one has 



e kin (/3, a,fi) - ^\a\ 2 p 



lim d |^W0,„, m (Va(W)) = dP(f3, a, p) 



The quantity between square brackets can be considered as the gauge invariant kinetic energy. 
Proof: Suppose that the theorem holds for a. = 0. We can then use a canonical transformation 
to obtain the result for arbitrary a as a consequence of (jlO.llfl : 

lo x (V(W)) = oo x ( la l (V(W))) = co~ x (V(W)) = 2e kin (P,p + ^cx 2 ) - dP((3,p+ ±\cx\ 2 ) 

By using (|10.12|) and (|1U.13|) . this is equivalent to the statement of the theorem. 

We now prove the theorem for a = 0. As the pressure is independent of the boundary conditions, 
we can use periodic boundary conditions to compute it, i.e., we choose A to be a d-dimensional 
torus. For t > 0, let tA be the torus rescaled by t. Then |tA| = t d \A\, and 

U t : L 2 (tA,dx) -» L 2 (A,dx) : (U t ip)(x) = t d/2 ^{tx) 

is unitary. The Laplacians on A and £Aare related as follows: 

u t A tA u; = t- 2 A A . 

This relation carries over to the kinetic energy in second quantization: 

U t H , tA u: = t~ 2 H , A . 
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where we have used the same notation for the corresponding unitary on the Fock space J-(L 2 (K)) 
with one-particle space L 2 (A). Similarly, one easily finds that the scaling behavior of the potential 
energy terms in the Hamiltonians is as follows: 

U t V tA (W)U* = V A ({W{t-)) . 

and for the particle number we have 

u t N tA u: = N A . 

By using these unitary equivalences we obtain 

P(/3, M ) = hm -LlogTr HL 2 (tA)) e-^ H ^+ v ^ w ^ N ^ 
A ^R d \tA\ v v " 



ififc W\ ' 0gTr ^Wje-*-^'"^-^" 



This shows that the last expression is independent of t. Setting its derivative in t = 1 equal to zero 
yields the following equation 

2e kin (/?,//) - lim -Lc^VaW) " dP{/3,n) = . 
A^R d |A| 



In order to close the dynamical equations, we need to express the expectation values of the 
currents iir,, given in (|3.5H3.8jl . in the states lo\ in terms of the expectations of the conserved 
quantities u 3 of (|2.6|) . 

Proposition 10.2 The expectations of the local currents w J k in a Gibbs state iv\ are given by 

^aC^J = aja k uj\(u° x ) + 5 k jP{\) 
^>.( w k,x) = a k uj x {ui) + a k P{\) 

where P(X) is the pressure defined in (|1U.1|) and u J x are the local densities of the five conserved 
quantities defined in (|2.6|) . With the definitions of (|2.15() and (|4.1|) . this is equivalent to A = w. 
Explicitly: q° = p, and 

^\{ w l,x) = a k p = q k 

^aO^J = a i a kP + Sj, kP = q j q k /q° + Sj, kP 
u x (wU = a k (q 4 + P)= q k (q 4 + P)/q° 
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Proof: The first equation, j = 0, follows directly from (|3.5|) . ()lU.7j) . and (|1U.11|) . The expressions 
for wi x ,j = 1, 2, 3, contain the virial of the potential W, which we can relate to the thermodynamic 
pressure by using the virial theorem, Theorem 1 1U.1> 



w Kx = VjOjVfcOa. - - W(x- y) 



2 L ' \x — y\ 

For j ^ k, the expectation of the second term vanishes as it changes sign under rotation over tt 
about the jth axis, which is a symmetry of the potential and the Gibbs states. Due to the rotation 
invariance of the potential, we also have W'(x)x/|x| = \7W(x). Therefore, the expectation of the 
second term in a Gibbs state u>\ is given by 

To treat the first term of wi , as well as the first two terms of w kx , we will transform these these 
terms to a frame where the Gibbs state has zero total moment, so that we can more easily use 
invarance under reflections in space. E.g., from (|l(J.llj) we get 

jd+V k a x ) = u)~ x {"f a {V ja+V k a x )) 

= uJ x (Vja x \7 k a x ) + aja k uj x (a x a x ) + iuo x (aja^V k a x - a k \7ja x a x ) 

As the total momentum has zero expectation in uj x , the last term vanishes for all j, k = 1,2, 3. By 
reflection symmetry and the defintion of the kinetic energy we have 

2 1 

u~ x {V jatV k a x ) = -e kin (/3,/i+ ^\a\ 2 ) 

By combining the above relations we obtain 

1, 



u\{^ k ,x) = aja k ujx(u x ) + -S jtk 



2e kin (/5,/i + ^H 2 ) - u\{y x ) 



= aja k uj x {u° x ) + -P((3,a,n) 

where, for the last equality, we have used the virial theorem and (|1U.13|) . 

To compute the energy current, wa(w£ J, we need to consider the following expectations: 

iujx(\ 7 kat a y a y a x ~ a x aya y V k a x ) 
iLUx(a x Vja,ya y a x - a x a+Vja y a x ) 
iux(\7 k a+Aa x - Aa+V k a x ) 
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Again, we use (|1U.11|) to relate these expectation to expectations in wr. The first expectation 
becomes: 

iuJ x (V k a^a^a y a x - a+a+a y Vka x ) + 2a k u~ x (a+a y v a y a x ) 
The first term of this expression vanishes by symmetry. In the same way we find 



iLU\(a x V ' jOy a y a x — a x a y V ja y a x ) = 2ajL0 x {a x a y a y a x ) 

We treat the third expression with similar arguments: 

iuj x (V k a^Aa x - Aa^V k a x ) = iu x (-f a ((V k a^ Aa x - Aa^V k a x )) 

= iu x (V k a^Aa x + 2iV fc a+(a • V)a x - \a\ 2 V k a^a x 

— ia k a x Aa x + 2a k a x (a ■ V)a x + ia k \a\ 2 a x a x ) 

+ complex conjugate 
= -a k u>~ x (2\a\ 2 a x a x + iV k a x V k a x - 2a x Aa x ) 

Then, by using integration by parts and reflection symmetry we get the following expression: 



(10.14) 
(10.15) 



iuj\{V k a^Aa x - Aa^V k a x ) = -Aa k 
Recall the expression for the energy current: 



\ a \ 2uJ l( u °x) + o e kin(A) 



(10.16) 



V k a^Aa x - AaJV k a x ^ + - j dyW(x - y) V k a^a y a y a x - a^a^a y V k a x 



W'(x - y) 



(x - y) k (x - y)j 



f - y\ 



y V J U I/ 1 



Using (|10.16j) . we see that the expectation of the first term in lu\ equals 



5-1 

ge kin (A) + -\ol\ 2 u x {u° x ) 



For middle term we use (|10.15j) and find 



a kU\(- I dyW(x - y)a x a y h a y a x ) 



Similarly, for the last term we get 



^a k uj x { j <h) 



W'{x-y) 



(x - y) k {x - y)j 



\x -y\ 



--^OL k uj\(v x ) 



P(A) - -e kin (A) 



where we have used the definition of v x (jl0.4fl and the virial theorem (Theorem 110. lj) . 

By combining the three terms and applying the relation (jlU.12|) one obtains the expression for 
oj\{w\ x ) given in the statement of this proposition. ■ 
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11 Appendix. The entropy inequality 

Our arguments rely in a crucial way on the following entropy inequality (|3.17j) : For any pair of 
density matrices 7 and uj, and for all self-adjoint h, and any 5 > 0, one has 

j(h) < ^logTr e 5h+losuJ + 5- 1 S(j\uj) (11.1) 

The inequality holds in the more general context of normal faithful states on a von Neumann 
algebra |17j . Here we give a proof for density matrices that emphasizes the connection with the 
variational principle of statistical mechanics. 

Proof: Let h be self-adjoint, and (3 > 0. The variational principle of statistical mechanics 18 
states that 

-^Tr e~ pH = inf [Tr jH - P^Sh)) 

p 7 

where the infimum is taken over density matrices 7, and S(j) := — Tr 7 log 7, is the von Neumann 
entropy of 7. For any non-singular density matrix uj, define H = — (/3 _1 (/i + logw), and take (3 = 5, 
use 

5(7 I u;) = Tr 7(log7 — logw) = — Tr 7log(w) — 5(7) , 
and rearrange the resulting inequality to obtain pi. 1(1 . | 



Equality in Ijll.lJ) holds if and only if 
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The inequality (|ll.lj) can also be turned around: 

5(7 I oj) < j(h) - logTr e^ 10 ^ , (11.2) 

and one can then take the sup over h to obtain a characterization of the relative entropy (as was 
done [T7]^: 

5(7 I u) = sup \j(h) - logTr e h+lo H (11.3) 

h L - 1 

with equality iff u> = e~ h /Tr e~ h , i.e., iff h = log + constant x 1. 

In contrast to the classical case, if log uj and h do not commute, we generally have 

logTr e h+loguJ + logTr we' 1 . 

However, due to the Golden-Thompson inequality, i.e., for any pair of self-adjoint A and B, 

Tre A+B < Tre A e B 

we still have 

Tr^h - log Tr uje h < 5(7 | uj) . 
Whenever uj and 7 do not commute, the equality will be strict for all h. 
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